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I INTRODUCTION  AND  SUMMARY 


The  Army  Materials  and  Mechanics  Research  Center  (AMMRC)  is  studying 
several  types  of  armor  penetration  problems  including  penetration  with 
rods  and  very  high  velocity  pellets.  A basic  element  in  the  study  is 
a computer  code  (HEMP*),  which  simulates  the  two-dimensional  axisymmetric 
wave  propagation  in  a penetration  experiment.  Present  day  two-dimensional 
Lagrangian  wave  propagation  computer  codes  such  as  HEMP  represent  the 
gross  plasticity  aspects  very  well,  and  yield  satisfactory  predictions 
for  cases  where  fracturing  of  the  armor  plate  is  not  a significant  part 
of  the  penetration  process.  The  predictions  are  less  satisfactory 
however,  for  cases  where  microcracking  and  shear  banding  dominate.  The 
latter  mechanisms  become  more  important  with  increasing  armor  plate 
hardness  and  increasing  projectile  hardness  and  velocity,  and  are 
enhanced  by  certain  projectile  geometries.  Thus  it  appeared  necessary 
to  account  for  fracture  in  a more  detailed  way. 

The  overall  objective  of  the  program  reported  here  was  to  impro/e 

the  predictive  capability  of  HEMP  for  penetration  calculations  by 

incorporating  into  it  recently  developed  dynamic  fracture  models,  the 

SRI  NAG  (Nucleation  And  Growth)  models.  These  models  for  brittle  and 
2-7 

ductile  fracture  provide  for  the  nucleation  of  cracks  or  voids, 
growth  of  these  flaws  during  the  period  of  tensile  loading,  and  the 
reduction  of  stress  and  strength  associated  with  the  developing  damage. 

The  models  appeared  to  be  appropriate  for  handling  the  penetration 
problem,  but  they  were  available  only  in  a one-dimensional  form.  There- 
fore, two  of  the  specific  objectives  of  the  study  were:  (1)  to  modify 
the  brittle  and  ductile  fracture  models  so  that  they  could  handle  two- 
dimensional  flow  and  (2)  to  incorporate  these  models  into  subroutines 
for  use  with  the  HEMP  code. 

To  be  able  to  use  a fracture  model  one  must  know  the  fracture 
parameters  appropriate  to  materials  of  interest.  Fracture  parameters 


9 


for  the  SRI  models  were  available  only  for  OFHC  copper,  1145  aluminum, 

Lexan  polycarbonate  (a  transparent  plastic),  and  several  grades  of 
beryllium.  Therefore,  a third  objective  of  the  study  was  (3)  to  conduct 
impact  experiments  in  an  armor  steel  and  derive  the  fracture  parameters. 
Included  in  this  objective  was  the  metallographic  examination  of  the 
targets  to  guide  in  derivation  of  the  fracture  models. 

During  this  program  the  above  objectives  were  largely  met.  The 
fracture  models  were  developed  and  incorporated  into  HEMP  in  several 
stages.  First  the  one-dimensional  models  were  simply  altered  to  account 
for  two-dimensional  flow.  For  brittle  fracture  this  required  constructing 
a new  description  of  the  crack  size  and  orientation  distributions.  An 
armor  steel,  XAR30,  was  selected  for  characterisation.  l'apered-f Iyer 
impacts  were  conducted  to  determine  the  compressive,  tensile,  and  fracture 
behavior.  Rod  impact  experiments  were  conducted  to  study  the  penetration 
behavior.  Posttest  metallographic  examination  was  made  of  the  sectioned 
targets  to  determine  qualitatively  the  nature  of  the  fracture  damage. 

For  the  tapered-f Iyer  impacts  the  cracks  were  also  analyzed  quantitatively 
to  determine  the  size  and  orientation  distribution  of  cracks.  From  these 
distributions  the  fracture  parameters  were  obtained  for  the  XAR30  armor 
steel,  (Table  V).  The  models  were  then  implemented  and  tested  in  a 
SRI  two-dimensional  code  similar  to  HEMP  and  used  to  simulate  successfully 
some  impact  problems  (Figure  24-27).  Finally,  the  models  were  modified 
to  account  for  higher  levels  of  damage.  The  models  were  then  delivered 
to  Mr.  John  Mescall  of  AMMRC  for  insertion  into  the  HEMP  code.  At  the 
time  of  writing  of  this  report,  the  models  for  ductile  and  brittle  fracture 
have  been  incorporated  into  HEMP,  and  test  calculations  have  been 
performed . 

The  remainder  of  this  report  is  organized  in  the  following  way. 

In  Section  II  the  fracture  model  in  two  space  dimensions  is  developed 
and  described  both  for  brittle  and  ductile  fracture  modes.  In  Section  III 
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the  experimental  work  is  described,  and  in  Section  IV  we  present 
computational  simulation  (with  the  developed  fracture  model)  of  three 
two-dimensional  impact  experiments.  Finally,  in  Section  V a summary  is 
given  of  the  results  of  the  program.  In  addition,  six  appendices  provide 
additional  details  concerning  the  computational  model  and  the  experiments 
performed. 


11 


II  MODELS  FOR  BRITTLE  AND  DUCTILE  FRACTURE 
IN  TWO  SPACE  DIMENSIONS 

A.  Introduction 

Computational  models  describing  brittle  and  ductile  fracture  were 
modified  to  include  two-dimensional  (planar  or  axisymmetric)  behavior 
and  were  applied  to  some  two-dimensionnl  wave  propagation  problems. 

In  these  computational  models  for  fracture,  damage  occurs  as  the 
nucleation  and  growth  of  small  voids  (ductile  fracture)  or  microcracks 
(brittle  fracture).  Nucleation  may  occur  physically  by  widening  of 
inherent  flaws  in  the  material,  cracking  of  hard  inclusions,  separating 
along  grain  boundaries,  or  by  other  mechanisms.  In  the  model,  however, 
nucleation  means  the  appearance  of  the  void  or  crack  of  an  observable 
and  easily  identifiable  size  on  photomicrographs  at  a scale  of  about  100X. 
This  nucleation  occurs  in  the  model  as  a function  of  stress  and  stress 
duration.  Following  nucleation,  the  voids  or  cracks  grow  at  a rate  that 
is  dependent  on  the  stress  level,  duration  of  loading,  and  the  size  of 
the  void  or  crack.  The  model  also  accounts  for  the  stress  reduction 
that  accompanies  the  development  of  damage. 

The  models  developed  in  this  work  were  for  incipient  damage  and  not 

for  full  separation.  The  brittle  fracture  model  has  been  extended  to 

8 

full  separation  on  a concurrent  project.  In  this  extended  model  the 
microcracks  coalesce  and  form  fragments.  The  ductile  fracture  model  had 
been  previously  verified  for  damage  resulting  in  10%  porosity  but  not 
for  greater  damage.  During  this  effort  the  program  was  modified  to  extend 
the  model  to  full  separation,  but  that  extension  has  not  yet  been  verified 
by  correlation  with  experiment. 

The  two  computational  models  of  fracture  are  implemented  in  subroutines 
that  may  readily  be  inserted  into  two-dimensional  Lagrangian  wave  propa- 
gation computer  codes.  While  the  material  is  undergoing  fracture,  these 
subroutines  are  called  instead  of  the  usual  equation-of-state  subroutines. 
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The  analytical  basis  for  the  two  models  is  presented  below.  These 

2-4 

analyses  are  somewhat  different  from  those  presented  earlier  because 
of  the  extension  to  higher  damage.  Details  of  the  computational  pro- 
cedures used  are  given  in  Appendix  A.  Implementation  of  the  models  into 
a computer  program  is  described  in  Appendix  B. 


B.  Ductile  Fracture 

The  ductile  fracture  model  was  formulated  on  the  basis  of  observations 

2 3 

of  ductile  fracture  in  soft  aluminum  and  copper.  These  observations, 

which  were  made  on  polished  cross  sections  of  targets  after  impact,  showed 

that  the  fracture  occurred  by  the  nucleation  and  growth  of  nearly  spherical 

voids.  The  observed  voids  were  measured  and  counted  and  assembled  into 

number-versus-radius  size  distributions.  These  surface  distributions  were 

2 

then  transformed  statistically  to  volumetric  distributions  with  the  BABS1 
computer  program.  A sample  set  of  void  size  distributions  is  shown  in 
Figure  1.  Curves  are  given  for  four  depths  within  the  semple;  the  maximum 
damage  is  at  the  plane  F01.  All  the  voliiretric  distributions  obtained 

..a**'” 

with  aluminum  and  copper  had  a form  that  could  be  approximated  by  the 
equation 


N (R)  = N exp  (-R/R  ) (1) 

go  1 

3 

where  N is  the  cumulative  number/cm  of  voids  with  radii  larger  than  R, 
g 

N is  the  total  number  of  voids  per  cubic  centimeter,  and  R^  is  a parameter 
of  the  distribution. 

The  total  void  volume  is  obtained  by  integrating  over  the  entire 
distribution. 


SP 


3 

= 8fiN  R 
o 1 


13 


CUMULATIVE  NUMBER  LARGER  THAN  R — numtoer/cm3 


FIGURE  1 OBSERVED  VOLUME  DISTRIBUTION  OF  VOIDS 
ON  THE  FREE  SURFACE  SIDE  OF  THE  SPALL 
PLANE  IN  AN  OFHC  COPPER  TARGET  AFTER 
A ONE-DIMENSIONAL  IMPACT:  EXPERIMENT  S24 


The  void  size  distribution  at  any  time  and  at  any  point  can  be  represented 

by  N and  either  R or  V . For  computational  purposes  N and  V are 
o 1 v o v 

selected . 

1.  Nucleation:  Nucleation  in  the  model  occurs  as  the  addition  of 
new  voids  to  the  existing  set.  These  new  voids  are  presumed  to  occur  in 

a range  of  sizes  with  a size  distribution  given  by  Eq.  (1).  At  nucleation, 

the  parameter  R equals  R , the  nucleation  size  parameter  (a  material 
1 n 

constant).  The  number  of  voids  nucleated  is  governed  by  a nucleation 
rate  function  that  was  derived  from  our  work  in  both  ductile  and  brittle 
materials . 

P - P 

N = N exp  ( -°)  P > P (3) 

o P no 

1 

= 0 P < P 

no 

where  N , P , and  P are  material  constants,  and  P is  the  tensile  pressure, 
o no  1 

The  constant  P is  the  threshold  for  nucleation. 
no 

The  void  volume  nucleated  in  a time  interval  *s  found  from 
Eqs.  (2)  and  (3)  to  be 

. 3 

AV  = 8tt  N At  R (4) 

n n 

2.  Growth:  In  the  model,  damage  increases  by  nucleation  of  new 
voids  and  by  growth  of  the  existing  voids.  In  our  studies  of  aluminum 
and  copper  it  was  found  that  growth  was  linearly  dependent  on  the  pressure 
level  and  current  void  size  so  that  the  growth  rate  R is 

P - P 

R = — £2  R • (5) 

4 11 

where  P is  the  tensile  pressure;  P is  the  threshold  pressure  for 

go 

growth,  and  1]  is  the  material  viscosity.  This  is  the  usual  form  for  a 
growth  law  in  a viscous  material  with  no  strength.  In  Reference  3 it 
is  shown  that  Eq.  (5)  is  accurate  for  small  voids,  but  for  larger  radii, 


- sSSjr 
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inertial  effects  reduce  the  velocity  below  that  given  by  Eq.  (5).  The 
growth  represented  by  Eq.  (5)  is  spherically  symmetric  because  the  void 
expands  equally  in  all  directions.  In  Reference  3 it  was  shown  that  Eq.  (5) 
is  an  appropriate  description  of  void  growth  in  material  with  strength 
undergoing  one-dimensional  planar  flow  as  well  as  for  spherical  flow.  The 
results  of  calculations  described  in  Appendix  C show  that  the  same  growth 
law  also  holds  for  conditions  of  high  shear  strain. 

The  growth  of  a void  during  a time  interval  At  is  obtained  by 
integrating  Eq.  (5)  to  obtain  the  new  radius  R, 


'P  - P 


R = R exp 
o 


go 


4T] 


At 


) 


(6) 


where  R is  the  radius  at  the  beginning  of  the  time  interval.  Since  every 


o 

void  in  the  distribution  grows  by  the  same  exponential  factor,  the  size 
parameter  R^  also  grows  according  to  Eq.  (6). 


R 


1 


(7) 


where  R is  the  size  parameter  at  the  beginning  of  the  time  interval. 
Then  the  new  void  volume  can  also  be  found  from  Eqs.  (2)  and  (7), 


3 

V = R 
v o 1 


(8) 


where  V = 877N  R , the  void  volume  at  the  beginning  of  the  time  interval, 
vo  o 10 

The  total  change  in  void  volume  is  the  sum  of  the  contributions 
associated  with  nucleation  and  growth.  Thus  the  total  volume  at  the  end 
of  the  interval  is 
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For  simplicity  in  the  computer  calculations,  a new  variable  is 
introduced  with  the  definition 


3 

T = — 

1 41) 


(10) 


3.  Pressure-Volume  Relation:  The  stress-strain  relations  for 

material  undergoing  fracture  account  for  the  presence  of  voids.  As 
usual,  the  stress  is  separated  into  pressure  and  deviatoric  components. 
The  deviatoric  stress-strain  relations  are  described  in  the  following 
subsection. 

The  pressure  is  related  to  the  specific  volume  and  internal  energy 
through  a combination  of  the  Mie-Grlineisen  equation  of  state  for  the 
solid  and  a relation  between  pressure  in  the  solid  and  average  pressure 
on  the  porous  material.  We  assume  that  an  average  pressure  in  the  solid 
material  can  be  computed  from  the  specific  volume  of  the  solid  and  the 
internal  energy  through  use  of  the  Mie-Grlineisen  equation: 


P = (Cp  + Dp2  + Sp3)  (1  - ~)  + Fp  E (11) 

s 2 s 

where  C,  D,  and  S are  constants,  T is  the  Griineisen  ratio,  E is  internal 

energy,  n is  the  solid  density,  and  p is  the  strain:  n = p /p  - 1,  where 
s so 

o is  the  initial  density  of  the  solid.  Here  we  neglect  all  nonlinear 
o 

terms  in  p because  the  solid  strains  and  stresses  are  all  small  during 
fracture,  and  obtain 


P = C / — — - l\  + rpE  (12) 

3 \ ?o  / 


The  pressure  computed  from  Eq.  (12)  is  necessarily  an  average  because  the 
actual  stress  states  will  vary  greatly  through  partially  fractured  material. 


The  average  pressure  on  the  gross  section  of  the-  fracturing  material 


can  now  be  related  to  the  pressure  in  the  solid  components  according  to 

9 

a relation  derived  by  Carroll  and  Holt  for  porous  material: 
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p = 


P V 
s s 


(13) 


where  P is  the  average  pressure  on  a section,  V is  the  specific  volume 

s 

of  the  solid,  and  V is  the  gross  specific  volume.  The  volume  V is  the 

sum  of  the  solid  volume  V and  the  volume  V associated  with  voids.  A 

s v 

combination  of  Eq?.(12)  and  (13)  serves  to  evaluate  the  pressure: 


\“o  Cs  } 


+ rpE 


(14) 


These  equations  may  also  be  inverted  to  obtain  expressions  for  the 

specific  volume  of  the  solid  in  terms  of  P or  P 

s 


V = 
s 


1 £E 
0q  + C 

P 

s 

1 + r 


(15) 


v - — — + I® 
s o cp  c 


(16) 


The  change  in  solid  volume  AV  is  related  approximately  to  the  change 

s 

in  solid  pressure  AP  by  differentiation  of  Eq. (15) 

s 


TAE  - V AP 
so  s 

AV  = V — 

s so  C/p  + te 


(17) 


where  V is  the  specific  volume  of  the  solid  at  the  previous  time  step 
so 

and  ae  is  the  change  in  internal  energy  during  the  time  step. 


4.  Deviator  Stress  Computation  in  Two  Dimensions:  The  deviator 
stresses  are  first  computed  elastically. 


q * — (7 ' + 2G 

1J  i jo 


(*ij  ‘ ~rZXi) 


(18) 


where  G is  shear  modulus,  Ae  is  the  strain  increment  tensor,  and 

A J 
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5a  is  the  Kronecker  delta.  In  the  fracture  routines  all  stresses 
(and  pressures)  are  positive  in  compression.  It  is  noted  here  that  in 
HEMP,  BFRACT , and  DFRACT, 


TXY 


and 


O 


12 


EXY 


c + f 
*12  *21 


2C 


12 


If  yielding  occurs,  then  the  deviator  stresses  are  computed  by  the 
visco-plastic  relation 


(19) 


where  d^  is  the  ij  component  of  the  plastic  strain  increment  and 
P ^ 

dy  is  the  scalar  "effective"  plastic  strain  increment.  Equation  (19) 

1 


is  solved  for  the  deviator  stresses  by  using  Wilkins'  procedure  to 
handle  the  first  term  on  the  right  side.  Then 


where 


= deviator  stress  from  Eq.  (18)  and 


(20) 


+ 2 


1 2 2 2 \ 
(ai2  + °23  + C13  ) 


The  damage  that  occurs  is  presumed  to  affect  both  the  yield  strength 

and  the  effective  shear  modulus  of  the  material.  The  modulus  is  reduced 

as  a function  of  the  developing  porosity  in  accordance  with  the  elastic 

10 

relations  of  Mackenzie.  His  formulation,  in  the  present  nomenclature, 
is 

G = G [1  - V PF]  (21) 

o L vK  J 
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where  G is  the  effective  shear  modulus,  G and  C are  shear  and  bulk 

o 

moduli  of  the  solid,  V is  the  specific  volume  of  voids,  p is  the  gross 
density,  and 


F 


3C  + 4G 

o 

9C  + 8G 


) 


o 


15 


(1  - y) 
(7  - 5y) 


(22) 


where  u is  Poisson's  ratio. 

The  factor  F varies  only  from  2.14  to  1.66  as  y ranges  from  o to  0.5. 

In  the  fracture  subroutine  F is  called  SMF  and  is  inserted  with  a value 

of  1.88  corresponding  to  V = 1/3. 

The  yield  strength  reduces  somewhat  more  rapidly  than  the  modulus 

as  the  porosity  increases.  Dynamic  calculations  of  void  growth  by  Seaman 
3 

et  al.  indicated  that  the  yield  strength  should  reduce  in  the  following 
way. 


Y = Y Tl  - 4V  pi 
oL  v 1 

This  expression  is  used  in  the  fracture  subroutine. 


(23) 


C.  Brittle  Fracture 
* 

The  brittle  fracture  model  was  formulated  on  the  basis  of  observations 

3,4  4 6 

of  fracture  in  Armco  iron,  beryllium,  novaculite  (a  fine-grained  quartz), 

5 

and  Lexan  polycarbonate  (a  transparent  plastic).  In  impact  experiments 

with  these  materials,  fracture  occurred  by  the  nucleation  and  growth  of 

microcracks.  These  cracks  were  measured  for  length  and  angular  orientation 

with  respect  to  the  direction  of  loading.  The  observed  cracks  were  then 

organized  into  groups  according  to  size  interval  and  angle  interval. 

These  surface  distributions  were  then  transformed  statistically  to 

3 

volumetric  distributions  in  size  and  angle  with  the  BABS2  computer 

program.  For  this  transformation  it  was  assumed  that  the  cracks  were 

penny-shaped  and  that  the  distribution  was  axisymmetric  around  the 

direction  of  propagation.  A sample  set  of  crack  distributions  is  shown 

* Fracture  was  termed  brittle  whenever  the  primary  damage  mode 

appeared  as  cracks.  In  many  cases,  such  as  Armco  iron,  the  cracks 
grew  in  what  is  often  termed  a "ductile"  manner. 
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in  Figure  2.  Here  the  angular  variation  has  been  suppressed,  so  the 
ordinate  is  the  total  number  of  cracks  larger  than  the  indicated  radius. 

The  volumetric  distributions  obtained  with  Armco  iron  had  the  same 
exponential  form  found  for  voids.  In  addition  to  a size  distribution, 
it  was  necessary  to  consider  an  orientation  distribution.  Furthermore, 
as  the  material  rotates  (as  it  does  at  the  rear  of  a plate  under  projectile 
impact)  the  cracks  rotate  with  the  material. 

Considering  the  requirements  for  size  and  orientation  distributions 
and  for  permitting  rotation,  the  fracture  model  was  constructed  with  an 
array  of  crack  orientations  or  "bins"  associated  with  each  computational 
cell.  Each  crack  bin  contains  penny-shaped  cracks  with  a specific 
orientation  and  with  a specific  size  distribution. 

The  size  distribution  for  each  bin  in  the  model  is  given  the 
following  analytical  form 

N1  = N1  exp  (-R/R1)  (24) 

go  1 

i th 

where  N is  the  total  number  of  cracks  per  cubic  centimeter  in  the  i 

i°  i 

bin,  N is  the  number  of  cracks  with  radii  greater  than  R,  and  R is  a 

B 1 th 

constant  giving  the  shape  of  rhe  crack  size  distribution  in  the  i bin. 

Each  bin  contains  cracks  normal  to  a specific  angle  cp  (in  the  x, 

y plane)  and  0 (in  planes  normal  to  the  x,y  plane)  as  shown  in  Figure  3. 

o o o o 

At  present  only  four  values  of  $ (0  , 45  , 90  , 135  ) and  two  of 
o o 

0 (0  , 90  ) are  used.  For  example,  bin  1 contains  all  cracks  with  cp 

o o o o 

between  -22.5  and  22.5  and  ip  between  45  and  135  . Note  that  cracks 

o o 

with  cp  between  157.5  and  202.5  are  included  in  this  bin.  Bin  5 contains 

o o 

all  cracks  with  0 less  than  45  or  greater  than  135  . In  the  model 

calculations  the  cracks  are  treated  as  though  they  all  had  orientations 

o o 

corresponding  to  the  center  of  the  bin  (cp  = 0 , 0 = 90  for  bin  1, 
o 

0=0  for  bin  5,  and  so  on.)  Experience  will  show  whether  five  is  a 
large  enough  number  of  bins. 
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FIGURE  2 CRACK  SIZE  DISTRIBUTION  IN  ZONES  NEAR  THE 
SPALL  PLANE  IN  AN  ARMCO  IRON  TARGET  AFTER 
A ONE-DIMENSIONAL  IMPACT:  EXPERIMENT  S25 


MA-2024-7 


FIGURE  3 CRACK  ORIENTATION  BINS  AND  DEFINITIONS  OF  0 AND  0 


In  two-dimensional  axisymmetric  or  planar  flow,  the  cells  may  rotate 

in  the  x-y  plane.  This  rotation  is  accounted  for  by  allowing  the  crack 

bins  to  rotate  with  the  material.  The  angular  rotation  is  given  by  the 

th 

variable  p . The  position  of  the  i bin  in  the  x-y  plane  is  then  given 
by 


rp.  = cp 


io 


+ P. 


(25) 


where  rp  is  the  current  angular  position  and  cp^  its  initial  position. 

At  the  end  of  a computation  p^  is  listed  so  that  the  bin  orientations 
can  be  related  to  the  fixed  x-y  grid.  The  crack  bin  rotation  is  independent 
of  the  initial  orientation  of  the  computational  cell  and  of  rotations  that 
the  cell  undergoes  before  fracture  begins. 

The  stresses  applied  normal  to  the  cracks  in  each  bin  are  determined 
by  the  usual  transformations  for  two-dimensional  problems.  The  stress  in 
the  Z direction  is  always  a principal  stress.  These  normal  stresses, 
which  govern  nucleation,  growth,  and  expansion  of  the  cracks  are 


a + (T 


a - 

x 


CP0 


COS  2m  + 


T si n 2m 
xy  ^i 


(for  0 = 90  ) 


(26) 


(for  0=0) 


(27) 


If  the  cell  rotates  by  an  angle  p^  (positive  counterclockwise),  then  the 

o 

angle  to  the  crack  group  becomes  cp  + p for  the  set  of  groups  at  0 = 90  . 

io  t 

Then  Eq.  (26)  is  used,  with  cp^  determined  by  Eq.  (25). 

The  cracks  are  presumed  to  open  elastically  to  the  value  given  by 
Sneddon^1 


_ 4(1  - u2) 
KE 


(28) 


where  6 is  one-half  the  maximum  separation  of  the  crack  faces  and  E is 
Young's  modulus.  The  crack  faces  form  an  ellipsoid  with  three  semi-axes 
6,  R,  and  R.  Then  the  volume  of  a crack  is 


24 


(29) 


v = 16(1  - U ) ,3^ 

lc  3E  cpifj 

The  volume  of  the  entire  crack  distribution  is  obtained  by  combining 
Eqs.  (24)  and  (29)  into  the  following  integral: 


• £•■■1 :/■ 

j j •'n 


drN1  exp  (-R/Ri)]dR 
L o 1 


lc 


dR 


32(1  - 


^ l < <»: 


,3  1 
1 


(30) 


1.  Crack  Nucleation:  Nucleation  in  the  model  occurs  as  the  addition 

of  new  cracks  to  the  existing  set.  These  new  cracks  are  presumed  to  occur 

in  a range  of  sizes  with  a size  distribution  like  Eq.  (24).  At  nucleation, 

the  parameter  R*  equals  R , the  nucleation  size  parameter  (a  material 
1 n 

constant).  The  number  of  cracks  nucleated  is  governed  by  a nucleation 
rate  function  similar  to  that  used  for  ductile  fracture: 

- a 


N = N exp 
o 


(31) 


where  N , a . and  a are  fracture  parameters  and  n is  stress  normal 
o no  1 

to  the  plane  of  the  cracks.  This  form  of  nucleation  function  resembles 

12 

the  relation  deduced  by  Zhurkov  for  the  rate  of  breakage  of  atomic  bonds. 

We  have  found  it  applicable  to  ductile  materials  and  also  to  such  diverse 

3 * 5 

brittle  materials  as  Armco  iron,  beryllium,  polycarbonate,  and  a three- 

13 

directionally  reinforced  quartz  phenolic  composite.  The  new  cracks  are 


* 4 

In  beryllium  it  was  found  that  the  deviator  stress  governs  nucleation 


and  not  the  stress  o' 
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nucleated  with  a range  of  sizes  such  that  the  number  greater  than  R is 


i 

aN  = N At  exp  (-R/R  ) (32) 

g o 

• i th 

where  N At  is  the  total  number  nucleated  in  the  i bin,  R is  the 

o 

nucleation  distribution  parameter,  and  At  is  the  time  step.  The  volume 
of  the  entire  nucleated  distribution  is  obtained  by  combining  Eqs.  (30) 
and  (32). 

. . i 


ve<=e  f 

i ^o 


dfN  /\t  exp  (-R/R  )]dr 

f 

lc  dr 


32(1  - V2)  At  R3 


(33) 


-E 
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If  the  material  is  isotropic  and  under  a uniform  tensile  stress, 

nucleation  of  cracks  occurs  with  equal  probability  in  any  direction. 

In  such  a case  the  number  of  cracks  assigned  to  each  bin  is  proportional 

to  the  solid  angle  subtended  by  the  bin.  The  fraction  of  the  total  solid 

angle  for  each  bin  is  called  FNUC  in  the  fracture  subroutine:  it  is  now 

set  for  isotropic  nucleation.  If  the  material  is  not  isotropic,  FNUC 

can  be  altered  to  reflect  the  observed  flaw  orientations. 

2.  Crack  Growth:  The  growth  law  derived  from  experimental  data 

4 

on  both  ductile  and  brittle  fracture  is: 


where  T^ 
Here  a 


“ = T (<y  - a >R 
dt  1 go 


(34) 


is  a growth  coefficient  and  (j  is  the  growth  threshold  stress, 

go 

is  treated  as  a constant  material  parameter,  but  in  some  cases 


go 

it  has  been  taken  as  the  critical  stress  for  crack  growth  according  to 
fracture  mechanics 
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* I_  V 

CT  = K 

ugo  4R  Ic 


(35) 


where  K is  the  fracture  toughness.  Since  a is  generally  very  small 
Ic  go 

for  impact  problems,  a in  Eq.  (34)  can  usually  be  taken  as  a small 

go 

constant. 

When  Eq.  (34)  is  integrated  over  a time  step  At  (holding  o constant), 

go 

the  final  value  of  the  radius  is 


R = R exp  [T  (a  - a )At] 

1 1 cpi/j  go  J 


(36) 


where  R is  the  radius  at  the  beginning  of  the  interval  and  o is  the 
1 cp^ 

average  stress  in  the  interval.  When  Eq.  (30)  is  combined  with  Eq.  (36) 
the  crack  volume  associated  with  growth  at  the  end  of  the  time  step  is 
found  to  be 


V1  = 32N1  (~-^)  (rV  a\  exp  [3T  (ff* , - ff  >At1 
g o E 1 1 cp0  g° 


(37) 


th 

The  total  number  of  cracks  in  the  i bin  at  the  end  of  the  time  step  is 

(38) 


i i • i 

N = N + N At 
1 o 


The  total  volume  V1  may  be  represented  as  the  sum  of  V*  and  V*  from 
t n g 

Eqs.  (33)  and  (37),  or  the  combination  of  cracks  may  be  described  by  a 

single  analytical  form  with  N*  cracks  and  a new  shape  parameter  R* 

1 2 


or 


i i i 

V = V + V 
t n g 


i i (1  - V ) i 3 i 

V = 32  N . - — (R  ) a , 

t 1 E 2 cpij) 


(39) 


(40) 
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Equating  the  two  expressions  for  V provides  a means  for  evaluating  R , 

t 2 

the  distribution  parameter  appropriate  to  the  end  of  the  time  step. 


i i 3 -i 

N (K  ) exp  [3T  (a  . 
,„i,3  o 1 1 c Dl!) 

(V  = p 

N 


. .i  3 
-a  ) it  + N /\t  R 
go  o 


(41) 


Now  the  damage  at  the  end  of  the  time  step  can  be  completely  characterized 

by  two  parameters,  N*  and  R* , obtained  from  Eqs.  '(38)  to  (41). 

1 2 

3.  Stress-Strain  Relations:  The  stress-strain  relations  for  material 
undergoing  brittle  fracture  are  the  same  as  those  for  ductile  fracture. 

The  void  volume  is  replaced  by  the  crack  opening  volume  in  the  equations. 
Because  the  crack  opening  volume  is  small,  the  usual  stress-strain  relations 
are  modified  very  little  by  the  damage. 

For  the  stress  calculation,  the  damage  is  treated  as  if  it  were 
isotropic.  Future  modifications  may  be  made  to  account  for  the  anistropic 
nature  of  the  damage. 

4.  Implementation  in  the  Computer  Program:  The  foregoing  analyses 
set  forth  the  basic  equations  describing  nucleation  and  growth  of  voids 
and  cracks,  and  the  damage-caused  modifications  of  the  stress-strain 
relations.  These  equations  are  used  directly  in  the  computer  subroutines. 
Further  derivations  required  for  the  numerical  procedures  in  the  subroutine 
are  described  in  Appendix  A.  The  subroutine  and  test  cases  for  it  are 
given  in  Appendix  B. 
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Ill  EXPERIMENTAL  INVESTIGATION 


The  purpose  of  the  experimental  program  was  to  examine  qualitatively 
and  quantitatively  the  nature  of  fracture  in  one  armor  steel.  The  steel 
selected  was  a high-hardness  armor  steel  designated  XAR30.  The  dynamic 
experiments  were  tapered-f Iyer  impacts  and  long  rod  penetration  tests. 

The  tapered-f Iyer  impacts  were  conducted  to  determine  both  the 
Hugoniot  (stress-strain  character  under  compression)  and  the  fracture 
(NAG)  parameters  in  tension.  These  impacts  were  instrumented  with  stress 
gages  at  the  rear  of  the  target  to  obtain  the  needed  Hugoniot  data.  The 
targets  were  sectioned,  and  the  crack  size  and  orientation  distributions 
were  measured;  from  these  distributions  the  NAG  fracture  parameters  were 
obtained. 

Long  rod  penetration  tests  were  conducted  to  provide  both  qualitative 
guidance  in  developing  the  two-dimensional  fracture  model  and  a test 
case  for  future  computer  simulation.  Besides  tests  on  XAR30  targets, 
penetration  tests  were  also  made  with  1145  aluminum  and  Armco  iron,  two 
materials  whose  fracture  behavior  was  already  well  characterized.  The 
aluminum  and  Armco  iron  tests  are  described  in  Appendix  D. 

For  further  insight  into  the  basic  mechanisms  associated  with 
penetration,  posttest  observations  were  made  on  targets  of  an  improved 
homogeneous  steel  (IHS)  following  penetration  experiments.  These  results 
are  described  in  Appendix  E.  The  micromechanics  of  the  fracture  behavior 
of  the  XAR30  and  IHS  steels  and  Armco  iron  are  compared  in  Appendix  F. 

A.  Materials 

W gn~hardness  armor  steel  made  by  Great  Lakes  Steel  Company  and 

United  States  Steel  Company  and  designated  as  XAR30  was  used  in  this 

14 

work.  The  chemical  compositions  reported  by  Hickey  are  given  in 

Table  I,  and  tensile  properties  are  given  in  Table  II.  The  microstructures 
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Table  II 


TENSILE  PROPERTIES  OF  XAR30  ARMOR  STEEL* 


0.2%  Yield 

Tensile 

Reduction 

Supplier  and 

Strength 

Strength 

Elongation 

in  Area 

Orientation 

(ksi) 

(ksi) 

(%) 

(%) 

Great  Lakes  Steel 


Longitudinal 

208 

262 

14.5 

49.8 

Long  transverse 

216 

253 

10.5 

39.8 

Short  transverse* 

211 

261 

10.0 

12.5 

U.S.  Steel 

Longitudinal 

201 

231 

11.5 

47.8 

Long  transverse 

204 

241 

11.0 

44.0 

* Short  transverse  properties  are  the  average  of  three  tests  performed 
at  SRI  and  described  more  fully  in  Appendix  E.  All  other  values  are 
the  average  of  two  tests  performed  by  Hickey. 
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t. 


of  both  materials  appear  to  be  a mixture  of  martensite,  bainite,  and 
retained  austenite.  Fracture  toughness  testing  by  Hickey  indicated 
that  the  United  States  Steel  material  was  slightly  tougher  than  the 
Great  Lakes  Steel  material. 

The  Vickers  hardness  of  the  material  was  measured  to  be  545  (corre- 
sponding to  R 52  and  BHN  509)  and  differed  by  less  than  3%  on  three 
c 

mutually  perpendicular  faces.  Small  unconventional  tensile  specimens 
were  prepared  from  the  half-in jh  plates  in  the  plate  thickness  direction, 
and  the  short  transverse  tensile  properties  were  measured.  The  results 
are  summarized  in  Table  II  and  described  in  more  detail  in  Appendix  E. 

The  strengths  and  elongations  ar>3  essentially  the  same  in  all  three 
directions  in  the  steel,  but  the  reduction  in  area  at  failure  is  much 
smaller  in  the  short  transverse  tests. 

B.  Flat  Plate,  Tapered-Flyer  Impact  Experiments 

Five  flat  plate  impact  experiments  were  performed  with  tne  gas  gun 
to  establish  the  dynamic  fracture  behavior  of  XAR30  armor  steel.  The 
experimental  configurations  are  given  in  Table  III.  The  objective  of 
the  first  two  experiments  was  to  measure  the  load  history;  the  objective 
of  the  remaining  three  experiments  was  to  produce  different  degrees  of 
fracture  damage  that  could  later  be  quantified.  As  discussed  in 
Section  IV,  the  measured  load  histories  are  correlated  with  the  fracture 
damage  to  obtain  parameters  that  describe  the  dynamic  fracture  behavior. 

The  experimental  arrangement  is  shown  in  Figure  4.  Specimen  disks 
1.5  inches  in  diameter  were  machined  to  various  thicknesses  (see  Table  III) 
and  mounted  in  the  target  plate.  The  edges  were  beveled  at  an  8 degree 
angle  to  allow  the  specimen  to  release  easily  from  the  target  plate 
upon  impact  and  fly  into  the  rags  of  the  catcher  tank. 

Tapered  flyer  plates  2\  inches  in  diameter  and  having  an  average 
thickness  half  that  of  the  specimen  were  machined  from  the  same  material 
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Table  III 


DYNAMIC  FRACTURE  EXPERIMENTS  ON  XAR30  ARMOR  STEEL 


Flyer 


Experiment 

Number 

Specimen 

Thickness 

(mm) 

Thickness 
at  Center 
(mm) 

Angle 
of  Taper 
(deg) 

Flyer 

Velocity 

(mm/usec) 

Remarks 

2024-1* 

7.62 

3.81 

5.7° 

0.451 

Full  spall; 
clear  gage 
record 

2024-2* 

10.1 

4.44 

5.7 

0.357 

Large  con- 
tinuous 
fracture; 
clear  gage 
record 

2024-3 

5.08 

2.24 

5.7° 

0.226 

Numerous 

micro- 

fractures 

2024-4 

12.1 

5.58 

6.84° 

0.200 

Numerous 
micro- 
f ractures 

2024-5 

2.54 

1.12 

2.3° 

0.259 

Numerous 

micro- 

fractures 


* Instrumented  with  back  surface  ytterbium  piezoresistive  stress 
transducers. 


and  mounted  on  the  front  of  a 6-inch-long  aluminum  projectile.  The 
angle  of  taper  varied  for  different  experiments  from  2.3  to  6.8  degrees 
as  shown  in  Table  III.  The  impacting  surfaces  v.ere  ground  flat  and 
parallel  to  within  0.0005  inch. 

The  projectile  was  accelerated  down  the  barrel  of  the  gas  gun  upon 
sudden  release  of  compressed  helium,  and  careful  alignment  of  the  flyer 
plate  and  specimen  resulted  in  flat  plate  impact.  The  impact  velocity 
was  measured  by  electrical  contacts  at  the  gun  muzzle. 

Since  the  stress  duration  in  the  specimen  varies  according  to  flyer 
thickness,  tapered  flyer  experiments  have  the  advantage  over  untapered 
flyer  experiments  in  that  a range  of  stress  durations  can  be  obtained  in 
one  experiment.  Since  the  extent  of  dynamic  fracture  damage  depends  on 
stress  duration  as  well  as  on  the  magnitude  of  the  stress,  a range  of 
damage  can  be  produced  in  a single  experiment. 

Analysis  of  the  stress  history  for  a tapered  flyer  experiment  is 
more  complicated  than  for  the  uniaxial  strain  conditions  resulting  from 
untapered  flyer  impact,  because  a shear  wave  arises  when  the  dilatational 
wave  reflects  from  the  inclined  rear  surface  of  the  flyer.  A two-dimen- 
sional (planar)  wave  propagation  code,  however,  was  used  successfully 
to  compute  the  stress  histories  and  damage.  This  analysis  is  described 
in  Section  IV. 

Wave  profiles  were  recorded  close  to  the  rear  surface  of  two  of  the 
specimens.  Ytterbium  piezoresistive  stress  transducers,  mounted  in 
3/8-inch  thick  blocks  of  C-7  epoxy  that  were  glued  to  the  specimens, 
produced  the  oscilloscope  traces  shown  in  Figure  5.  Both  records  show 
the  Hugonlot  elastic  limit  (HEL),  a flat-topped  loading  wave,  and  a 
clear  fracture  signal.  The  HEL  mepsures  the  yield  strength  of  the 
material  under  dynamic  uniaxial  strain  conditions  and  is  useful  in 
specifying  the  constitutive  equation  of  the  armor.  The  flat-topped 
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loading  wave  indicates  that  impact  planarity  was  good  and  that  signi- 
ficant attenuation  did  not  occur  before  fracture.  The  second  peak  in 
the  records,  known  as  the  fracture  signal,  is  caused  by  a reloading 
recompression  of  the  transducers  by  waves  emanating  from  the  internal 
fracture  surfaces  as  they  form  and  grow. 

Shock  absorbing  materials  were  located  in  the  path  of  the  impacted 
specimen  to  decelerate  it  gradually  and  prevent  subsequent  shock  loads. 

The  recovered  specimens  were  then  sectioned  to  reveal  the  internal 
fracture  damage  as  shown  in  Figure  6.  Specimens  2024-3  and  2024-4  were 
sectioned  so  as  to  produce  a cross  section  in  the  direction  of  maximum 
taper  (see  insets  in  Figures  7 through  11;  Specimen  2024-5  was  sectioned 
parallel  to  the  direction  of  taper  at  three  locations  (see  inset  in 
Figures  12  and  13). 

The  fracture  damage  in  Specimens  2024-3,  2024-4,  ai  2024-5  was  analyzed 
quantitatively  by  counting  and  measuring  the  traces  of  the  microfractures 
on  the  polished  section  surfaces.  Measurements  were  made  by  using 
a large  area  record  reader  (LARR)  with  which  an  operator  positions  a 
cross  hair  on  one  end  of  a crack  trace,  pushes  a button  to  record  the 
coordinates,  and  repeats  the  process  for  the  other  end.  A simple 
computer  program  uses  these  data  to  compute  the  length,  orientation, 
and  position  within  the  specimen  of  the  trace.  If  desired,  a similar 
procedure  can  be  used  to  obtain  crack  widths.  The  size  distributions 
on  the  sectioned  surfaces  are  then  converted  to  actual  crack  size 
distributions  per  unit  volume  by  means  of  a statistical  transformation 
implemented  in  the  BABS2  computer  code.  This  procedure  is  described  in 
detail  in  Reference  3. 

The  results  of  this  quantitative  damage  analysis  are  presented  in 
Figure  7 through  13,  which  show  crack  size  distribution  curves  for 
various  positions  on  a cross  section  for  several  cross  sections  in 
Specimens  2024-3,  2024-4,  and  2024-5,  No  fracture  damage  was  observed 
for  section  C-J  of  Specimen  2024-5. 
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IMPACTI 


(■I  SPECIMEN  2024-4,  SECTION  A-B 
(Surface  ii  parallel  to  taper  of  flyer) 


lb)  SPECIMEN  2024-5,  SECTION  A A 

(Surface  b perpendicular  to  taper  of  flyer) 
MP-2024-1 6 


FIGURE  6 POLISHED  CROSS  SECTIONS  SHOWING  INTERNAL  SHOCK-INDUCED  FRACTURES 


CUMULATIVE  CRACK  CONCENTRATION  — nuirtar/cm' 


0.02  0.04  0.06  0.08 

CRACK  RADIUS  — cm 


MA-2034-11 


FIGURE  9 CRACK  SIZE  DISTRIBUTIONS  ON  SECTION  A-B,  FILE  1 OF  SPECIMEN  2024-4 
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CUMULATIVE  CRACK  CONCENTRATION  — numbar/cm 


CUMULATIVE  CRACK  CONCENTRATION  — ngmNr/cm 


C . Rod  Impact  Experiments 


For  the  rod  impact  experiments,  square  specimens  10  or  15  cm  on  a 
side  were  cut  from  the  half-inch  armor  steel,  and  the  front  and  back 
surfaces  were  ground  smooth  and  parallel.  The  projectiles  were  cylinders 
1.03  cm  in  diameter  by  2.03  cm  long  made  of  drill  rod  heat  treated  to 
various  hardness  levels  and  fitted  into  1.16-cm-diameter  Lexan  poly- 
carbonate sabots. 

The  experiments  were  conducted  with  a remotely  fired,  propellant- 
activated  gun,  Figure  14.  The  gun  consists  of  a 1917  Enfield  action 
fitted  with  a heavy  barrel,  which  is  chambered  for  the  1.16  Winchester 
magnum  cartridge  and  smooth-bored  to  1.16-cm  diameter.  Strain  gages 
were  mounted  at  three-  locations  on  the  barrel — over  the  chamber,  at  30 
cm,  and  at  38  cm  from  the  breech — and  were  routinely  monitored  during 
every  experiment.  The  pressure-versus-time  data  were  used  in  designing 
powder  charges  for  desired  projectile  velocities  and  in  ensuring  that 
safe  pressures  were  not  exceeded  in  the  gun. 

Projectile  velocity  was  determined  either  from  a calibration  curve 
relating  velocity-at-impact  to  powder  charge  or  by  measuring  the  times 
between  successive  cutoffs  of  three  light  beams  by  the  projectile  as  it 
emerged  from  the  muzzle.  Fiber  optics,  mounted  in  an  aluminum  muzzle 
extension,  were  used  for  transmission  and  pickup  of  the  light  beams. 

The  muzzle  extension  also  served  as  a projectile  guide,  ensuring  normal 
impact  of  the  projectile  on  the  target.  Since  over  half  of  the  projectile 
length  is  still  in  the  muzzle  extension  at  the  moment  of  impact,  the 
muzzle  extension  is  slotted  to  relieve  the  gas  pressure. 

Details  of  the  five  rod  impact  experiments  are  given  in  Table  IV,  and 
views  of  polished  cross  sections  of  the  specimens  through  the  point  of 
impact  are  shown  in  Figure  15.  No  fracture  damage  and  very  little 
deformation  occurred  at  0.64  mm/usec.  At  higher  velocities  a large 
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FIGURE  14  EXPERIMENTAL  FACILITY  FOR  LONG-ROD-IMPACT  TESTS 


FIGURE  15  POLISHED  AND  ETCHED  SECTIONS  THROUGH  8.66-mm-THICK  XAR30 
ROLLED  HOMOGENEOUS  STEEL  PLATES  SHOWING  THE  EFFECT  OF 
THE  VELOCITY  OF  THE  IMPACTING  ROD  ON  THE  DAMAGE  PATTERN 
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Table  IV 


ROD  IMPACT  EXPERIMENTS  ON  XAR30  ARMOR  STEEL 


Experiment 

Number 

Specimen 

Thickness 

(mm) 

Proj ectile 

Velocity 

(mm/usec) 

Remarks 

31 

8.66 

1.0 

Penetration;  gross  back 
surface  scab. 

32 

8.66 

0.85 

Penetration;  large  mid- 
plane crack. 

33 

8.66 

0.72 

Penetration. 

34 

8.66 

0.64* 

No  penetration. 

35 

8.66 

0.67* 

Penetration. 

* Optically  measured. 
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(a)  FULLY  SCABBED  BACK  SURFACE 


(b)  BACK  SURFACE  SCAB 


FIGURE  16  APPEARANCE  OF  XAR30  TARGET 
IN  LONG  ROD  EXPERIMENT  31 
AFTER  IMPACT  AT  1.0  mm/jusec 

50 


central  crack  formed,  Figure  15(c)  and  (d),  which  at  1.0  mm/usec  caused 


a large  scab  of  the  armor  steel  to  break  free,  Figure  15(e)  and  Figure 


16.  Figure  16  shows  the  appearance  of  the  scabbed  back  surface,  the 


scab,  and  the  plug  from  Experiment  31. 


Penetration  in  these  experiments  is  classified  as  occurring  in  a 


plugging  mode.  Close  examination  of  the  surfaces  shown  in  Figure  15 


shows  the  existence  of  narrow  bands  that  appear  white  when  etched  in 


nital  or  Vilella’s  reagent  (Figure  17).  These  bands  are  known  as 


adiabatic  shear  bands  and  act  as  preferred  cracking  paths.  Thus  the 


penetration  process  in  XAR30  armor  steel  includes  the  formation  of 


highly  localized  regions  of  intense  shear,  which  fail  and  result  in  the 


liberation  of  the  plug-like  segment  of  material  from  the  armor  plate. 


The  plugs  from  these  experiments  were  also  sectioned,  polished, 


etched,  and  examined  with  a microscope.  A small  meniscus-shaped  region 


of  different  etching  response,  Figure  18,  was  observed  directly  beneath 


the  impact  surfaces  of  Specimen  31,  indicating  that  the  0'bcc  ~ €hcp 


polymorphic  phase  transformation  had  taken  place.  This  zone  was  not 


observed  in  the  other  specimens.  It  was  recently  shown  that  the 


lower  boundary  of  this  meniscus  corresponds  to  an  isobar  of  about  130 


kbar.  It  was  also  found  that  the  occurrence  of  this  phase  change  signi- 


ficantly alters  the  stress  history  and  hence  the  fracture  damage  in  the 


target  material. 


Specimen  31  also  exhibited  several  cracks  parallel  to  the  impact 


surface  at  about  mid-thickness,  Figure  18.  Only  one  such  crack  was 


* Adiabatic  shear  bands  are  narrow  regions  of  highly  localized  large 
plastic  shear  strains.  The  heating  accompanying  the  shear  deformation 
may  increase  the  temperature  high  enough  to  cause  solid  phase  trans- 
formations or  even  melting.  Very  rapid  quenching  of  this  material 
follows  because  the  large  volume  of  adjacent  material  in  intimate 
contact  with  the  narrow  band  conducts  the  heat  away  at  high  rates.  In 
the  present  instance  a transformation  to  austenite  followed  by  another 
rapid  transformation  to  martensite  probably  occurred.  The  white  etching 
response  in  5%  nital  is  consistent  with  the  existence  of  martensite. 
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i 


(bl  SPECIMEN  35 


MP-2024-27 

FIGURE  17  ENLARGED  VIEWS  OF  SECTIONS  THROUGH 
SPECIMENS  35  AND  32  NEAR  THE  ZONE 

OF  PENETRATION  SHOWING  THE  ADIABATIC  * 

SHEAR  BANDS 
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FIGURE  18  ENLARGED  CROSS  SECTION  OF  THE  PLUG  FROM  THE  XAR30  TARGET 
IMPACTED  IN  LONG  ROD  EXPERIMENT  33 


i 

i 


< 
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observed  in  Specimen  32;  none  was  observed  in  the  other  specimens. 
Adiabatic  shear  bands  bordered  the  edges  of  all  plugs.  A single  shear 
band  running  inward  from  the  sides  at  mid -thickness  was  observed  in 
Specimens  31  and  32. 
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IV  APPLICATIONS  OF  THE  DUCTILE  AND  BRITTLE  FRACTURE  MODELS 


The  two-dimensional  ductile  and  brittle  fracture  models  were 
applied  to  the  simulation  of  damage  in  several  metals.  The  brittle 
fracture  parameters  for  XAR30,  an  armor  steel,  were  derived  from  impact 
experiments,  as  described  in  the  previous  section.  Then  model  calcula- 
tions were  performed  to  simulate  two-dimensional  damage  in  a target  of 
the  material.  To  further  demonstrate  the  models,  simulation  calculations 
were  made  for  tapered-f Iyer  impacts  in  1145  aluminum  (ductile  fracture) 
and  in  Armco  iron  (brittle  fracture). 

The  first  step  in  applying  a fracture  model  to  a material  is  to 

determine  the  fracture  parameters  governing  nucleation  and  growth  of 

damage.  The  fracture  parameters  for  XAR30  armor  steel  were  computed 

from  the  observed  damage  in  three  tapered-f Iyer  experiments.  The 

following  procedure,  which  was  used  to  calculate  the  parameters,  was 

4 

similar  to  that  developed  for  beryllium.  The  average  crack  nucleation 

rate  (total  number  of  cracks  divided  by  nominal  duration  of  the  tensile 

stress)  was  plotted  versus  peak  stress  in  tension  to  find  the  nucleation 

• 

threshold  stress  (cr  ) and  the  other  nucleation  parameters  (ff  and  N ). 

no  1 o 

c 

The  shape  parameter  R^  of  the  observed  distribution  was  plotted 

versus  tensile  impulse  (peak  tension  times  duration  of  the  tension)  to 

c 

determine  the  nucleation  size  (R  ) and  the  growth  rate  constant  (T  ). 

o 1 

After  the  initial  estimates  of  all  five  parameters  were  determined 
from  plots,  trial  one-dimensional  calculations  were  performed  to  approxi- 
mate the  impact  conditions  at  several  points  in  the  target.  These 
calculations  were  repeated  with  different  fracture  parameters  until  the 
computed  and  measured  damage  compared  satisfactorily.  Then  a two- 
dimensional  calculation  was  performed  to  simulate  the  entire  impact.  It 
was  not  necessary  to  modify  the  parameters  further  and  repeat  the  two- 
limensional  simulation.  The  fracture  parameters  found  for  XAR30  armor 
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are  listed  in  Table  V.  The  parameters  for  Armco  iron  and  aluminum, 

2-4 

obtained  on  earlier  projects,  are  also  included  in  the  table. 

Some  of  the  results  of  the  simulations  for  three  tapered-f Iyer 
impacts  are  shown  in  Figures  19  through  21.  These  figures  show  the 
crack-size  distributions  on  the  planes  of  maximum  damage  at  each  section. 
The  crack  size  distributions  are  exponential  in  the  computations  and 
therefore  appear  as  straight  lines  in  these  figures.  The  computed 
distributions  are  high  in  some  case  and  low  in  others  but  the  agreement 
is  considered  satisfactory. 

Two  fracture  calculations  were  performed  with  the  ductile  and 
brittle  fracture  models  to  examine  their  capabilities.  Each  was  the 
simulation  of  a tapered-flyer  experiment  in  a well-characterized 
material:  1145  aluminum  (ductile)  and  Armco  iron  (brittle).  These 

tapered-flyer  experiments  were  selected  instead  of  projectile  impacts 
because  no  large  distortions  of  the  Lagrangian  mesh  occur  during  a 
calculation.  Thus  it  was  possible  to  test  the  fracture  models  without 
the  difficulties  involved  in  rezoning  and  construction  of  slide  lines 
that  would  be  required  for  more  complex  geometries. 

The  configuration  for  the  tapered  flyer  experiments  was  shown 
schematically  in  Figure  4,  and  dimensions  for  the  sample  cases  are 
shown  in  Figure  22.  For  the  aluminum  impact  calculations,  the  width 
of  the  flyer  and  target  were  foreshortened  to  5,328  mm  to  minimize  the 
number  of  computational  cells  required.  In  the  computations  the 
lateral  boundaries  were  allowed  no  horizontal  motion,  thus  approximately 
simulating  the  conditions  in  the  central  region  of  a tapered-flyer 
target. 

The  first  wave  traveling  through  the  target  after  the  impact  was 
a compressive  wave.  This  wave  was  reflected  from  the  free  rear  surface 
of  the  target  as  a zero-stress  rarefaction  (waves  4 in  Figure  23).  This 
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CRACK  RADIUS  — cm 


MA-2024-17 

FIGURE  19  COMPARISON  OF  MEASURED  AND  COMPUTED  DAMAGE  ON 
THE  PLANES  OF  MAXIMUM  DAMAGE  IN  TAPERED-FLYER 
IMPACT  EXPERIMENT  2024-3 


CUMULATIVE  CRACK  CONCENTRATION  — number/cm' 


FIGURE  21  COMPARISON  OF  MEASURED  AND  COMPUTED  DAMAGE  ON 

THE  PLANES  OF  MAXIMUM  DAMAGE  IN  TAPERED-FLYER  4( 

IMPACT  EXPERIMENT  2024-5 
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(b)  ARMCO  IRON  FLYER  AND  TARGET,  SHOT  SI 
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FIGURE  22  CONFIGURATIONS  FOR  TAPERED-FLYER  IMPACT  EXPERIMENTS 


61 


X ti)in  ll&iv 


POSITION,  x 


MA-2161-40 

FIGURE  23  DISTANCE  TIME  PLOT 
SHOWING  WAVE  PATHS 
AND  COMPRESSIVE  (C) 
AND  TENSILE  (T) 
REGIONS  IN  A ONE- 
DIMENSIONAL IMPACT 


rarefaction  interacted  with  the  rarefaction  from  the  free  surface  of 
the  flyer  (waves  3)  to  produce  tension  in  the  target,  as  shown  by  the 
region  labelled  T in  Figure  23.  The  vertical  length  of  the  tensile 
region  is  the  duration  of  the  tension  at  a point  in  the  target.  This 
duration  is  a linear  function  of  the  flyer  thickness.  Therefore  in  the 
tapered-f Iyer  experiment,  the  peak  tensile  stress  is  constant,  but  the 
duration  varies  in  the  direction  of  the  taper. 

These  simulations  are  given  as  sample  problems  in  Appendix  B. 

The  computed  results  are  depicted  here  as  three-dimensional  surfaces 
in  which  the  damage  functions  are  the  amplitudes  above  the  plane  of  the 
sections  through  the  targets.  The  damage  in  the  aluminum  target  is 
shown  in  Figures  24  and  25.  Figure  24  shows  that  the  line  of  maximum 
void  volume  is  nearly  at  the  middle  of  the  target  and  that  the  amount 
of  void  volume  increases  toward  the  high  damage  end.  The  high  damage 
end  is  the  part  of  the  target  struck  by  the  thicker  portion  of  the 
flyer.  Similar  results  are  observed  with  the  crack  concentrations 
plotted  in  Figure  25.  We  note  the  damage  functions  present  a fairly 
smooth  surface,  as  is  expected  from  a calculation  if  there  are  no 
stability  problems. 

The  crack  volume  and  crack  concentrations  computed  for  the  Armco 
iron  impact  are  shown  in  Figures  26  and  27.  Note  that  both  the  amount 
of  damage  and  the  extent  through  the  thickness  increase  towards  the 
high  damage  end  of  the  target. 

Quantitative  analyses  of  the  experimental  damage  was  not  performed 
for  these  two  tapered-f Iyer  impact  experiments.  However,  a qualitative 
comparison  can  be  made  between  the  damage  observed  and  that  computed. 

The  aluminum  Impact  calculation  compares  well  with  the  observed  ampli- 
tude of  damage  and  the  width  of  the  damaged  region.  The  Armco  iron 
calculation  indicates  about  the  right  level  of  damage,  but  the  width 
of  the  damage  region  appears  broader  than  that  seen  in  the  experiments. 


63 


3.25  0 

DEPTH  IN  TARGET  — cm 


MA-2024-21 


FIGURE  24  COMPUTED  VOID  VOLUME  THROUGHOUT  THE  1145  ALUMINUM  TARGET 
AFTER  TAPERED-FLYER  IMPACT  EXPERIMENT  S4 


FIGURE  25  COMPUTED  NUMBER  OF  VOIDS  THROUGHOUT  THE  1145  ALUMINUM  TARGET 
i AFTER  TAPERED-FLYER  IMPACT  EXPERIMENT  S4 


FIGURE  26  COMPUTED  CRACK  VOLUME  THROUGHOUT  THE  ARMCO  IRON  TARGET 
IN  TAPERED-FLYER  EXPERIMENT  SI  AT  2.069  fmc 
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CRACKS  — number /cm 


This  lack  of  agreement  suggests  a need  to  allow  nonelastic  crack 
opening  relations  for  the  iron  to  permit  damage  to  reduce  the  stress 
more  rapidly.  Such  inelastic  opening  certainly  occurs  because  the 
cracks  remain  open  after  the  impact.  However,  appropriate  opening 
relations  have  not  yet  been  developed. 

In  summary,  the  brittle  fracture  parameters  for  XAR30  armor  steel 
were  derived  from  three  tapered-f Iyer  impact  experiments.  The  ductile 
and  brittle  fracture  models  were  then  used  to  perform  sample  two- 
dimensional  impact  calculations  to  simulate  damage  in  1145  aluminum  and 
Armco  iron.  The  results  of  the  computations  in  the  three  materials 
show  that  the  models  can  simulate  satisfactorily  the  ductile  and  brittle 
damage  in  two-dimensional  problems. 


Appendix  A 

DERIVATION  OF  EQUATIONS  FOR  FRACTURE  SUBROUTINES 

The  detailed  equations  required  to  implement  the  ductile  and  brittle 
fracture  models  in  wave  propagation  computer  programs  are  derived  in 
this  appendix.  These  models  are  the  bases  of  the  subroutines  DFRACr 
and  BFRACT,  which  replace  the  usual  equation-of-state  subroutines  when 
ductile  or  brittle  fracture  begins.  Thus,  during  the  early  part  of  an 
impact  calculation  while  the  material  is  in  compression,  the  usual 
equation-of-state  subroutine  is  called  to  compute  the  stress  for  each 
cell  at  each  time  increment.  Later,  when  the  computed  tensile  stress 
exceeds  the  nucleation  threshold,  DFRACT  or  BFRACT  is  called.  These 
routines  compute  the  pressure  and  stress  from  given  values  of  the  strain 
increments,  internal  energy,  and  other  parameters.  In  addition,  these 
subroutines  nucleate  cracks  or  voids  and  permit  existing  cracks  or  voids 
to  grow.  Once  a fracturing  routine  is  called  to  compute  stress  for  a 
computational  cell,  the  usual  equation-of-state  subroutine  is  never 
called  again  for  that  cell.  Thus  the  presence  of  damage  is  accounted 
for  in  subsequent  recompressions  or  extensions  by  these  fracture  sub- 
routines. In  the  model  the  recompression  wave  does  not  compact  voids 
or  reduce  the  number  of  voids  or  cracks;  it  simply  compresses  the  solid 
material.  Thus  the  model  treat3  correctly  only  a low  amplitude  re- 
compression wave.  When  a later  tension  occurs,  the  damage  again  begins 
to  increase  by  nucleation  and  growth. 

The  stress  and  damage  quantities  are  related  by  strongly  nonlinear 
equations  so  that  an  iteration  solution  procedure  is  required.  The 
development  of  this  iteration  procedure  and  of  equations  for  the  stress 
estimates  used  to  start  the  iterations  is  described  in  this  appendix. 
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Iteration  Procedure 


The  subroutine  is  organized  around  an  iteration  procedure  to 
determine  simultaneously  the  damage  quantities  and  the  stresses.  To 
minimize  the  number  of  iterations,  a complex  estimating  procedure  is 
derived  for  determining  the  starting  value  for  each  iteration.  Methods 
for  making  the  required  estimates  are  outlined  in  the  second  subsequent 
section  of  this  appendix. 


The  subroutines  are  provided  with  an  internal  energy  E and 
strains  Af  , as  well  as  the  values  of  E,  o,  pressure,  and  stresses  at 
the  previous  time  step,  and  are  asked  to  provide  pressure  and  deviator 
stresses  at  the  current  time.  These  new  stress  quantities  are  a function 
of  the  changes  in  energy  and  density,  and  also  of  the  growing  damage. 

Thus  we  require  the  simultaneous  solution  of  the  following  system  of 
equations: 

Damage  = f , (P,  c r'  , or'  , T , <r'  ) 

1 xx  yy  xy  00 


O = f (Damage,  AE,  Ap) 


(Al) 


where  P = pressure 


a'  , a'  , o'  = deviator  stress 
xx  yy  00 


T = shear  stress  on  the  x-y  plane 
xy 

(7  = any  stress 


The  iteration  process  for  solving  the  system  of  Eq.  (Al)  requires 
(1)  an  estimate,  (2)  a computation  of  all  quantities  including  one  with 
which  the  accuracy  of  the  estimate  can  be  tested,  and  (3)  a test  for 
convergence.  The  change  in  solid  volume  AV^  was  chosen  as  the  initial 
parameter  for  beginning  an  iteration  because  all  the  quantities  can  be 
computed  from  that  one  estimate.  Convergence  is  based  on  AV,  the  overall 
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volume  change.  The  iteration  process  contains  the  following  six  steps, 
which  provide  the  framework  of  the  subroutine. 

(1)  Estimate  AV  . 

s 

(2)  Compute  the  pressure  and  deviator  stresses. 

(3)  Compute  the  crack  or  void  volume,  V , from  the  growth, 

va 

nucleation,  and  expansion  laws. 

(4)  Compute  the  total  volume  change  AV  from  the  change  in  crack 

a 

or  void  volume  AV  and  the  change  in  solid  volume  AV  . 

va  s 

(5)  Compare  AV  and  AV  and  terminate  the  iterations  if  the 

a 

comparison  is  satisfactory. 

(6)  Reestimate  AV  and  return  to  step  (2). 

s 

The  accuracy  requirements  for  determining  convergence  were  developed 
from  trial  computations  with  Armco  iron  and  a program  to  test  BFRACT. 
Several  density  step  sizes  and  accuracy  controls  were  used.  It  was 
found  that  the  results  can  be  made  independent  of  step  sizes  to  a 
precision  on  stress  of  0.01  kbar: 

(1)  If  step  size  in  density  is  such  that  the  maximum  change  in 

stress  could  be  0.33  rr  , that  is,  one-third  the  nucleation 

no 

threshold  stress. 

5 

(2)  If  the  ratio  of  (AV  - AV  )/V  is  less  than  2 x 10  . 

a s 

If  the  iterations  do  not  converge  in  ten  tries,  an  abort  procedure 
is  provided.  Normally  the  convergence  improves  as  the  density  step 
size  decreases.  Therefore,  in  the  abort  procedure,  the  step  size  is 
decreased  and  the  calculations  are  repeated.  If  convergence  is  still 
not  achieved,  a message  is  printed  and  the  calculations  continue. 


71 


For  compressive  stresses  there  is  no  crack  volume  and  void  volume 

is  held  constant:  then  the  initial  estimate  of  AV  is  exact  and  no 

s 

iteration  is  required.  For  small  amounts  of  damage,  the  first  estimate 
is  usually  accurate  enough  to  provide  convergence  on  the  first  iteration 
cycle.  However,  a3  the  increase  in  crack  or  void  volume  becomes  comparable 
to  the  imposed  total  volume  changes,  the  number  of  iterations  increases. 
Even  at  large  damage,  convergence  usually  occurs  in  three  to  five 
iterations. 

Pressure  Estimation 

For  beginning  each  iteration  it  is  necessary  either  to  estimate 
the  specific  volume  directly  or  to  estimate  the  pressure  and  derive 
the  specific  volume  from  the  pressure.  The  initial  estimate  of  pressure 
is  made  as  accurately  as  possible  to  minimize  the  number  of  iterations 
required.  Estimates  are  made  on  the  assumption  that  the  pressure  is 
determined  entirely  by  strain  and  changes  in  internal  energy,  by 
nucleation  of  new  cracks  or  voids,  or  by  a combination  of  strain, 
expansion,  and  growth  of  cracks  or  voids. 

For  strain  and  changes  in  internal  energy  only,  Eq.  12  (See  Section 

II  of  this  report)  is  used,  with  p computed  by  assuming  that  all 

r s 

strain  is  taken  in  changing  the  solid  specific  volume,  i.e.,  AV  = AV. 

s 

The  estimate  of  pressure  in  the  solid  is 


where  V 

so 

time  step. 


P = C 
se 


(o  (V  + AV)  " 1) 
\ o so  / 


r e 


V + AV 
so 


(A2) 


is  the  specific  volume  of  the  solid  at  the  beginning  of  the 


If  there  are  few  cracks  or  voids  in  the  cell,  the  initial  pressure 
step  beyond  the  nucleation  threshold  may  be  governed  by  nucleation. 

For  cracks  the  pressure  estimate  based  on  nucleation  is  obtained  by 
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setting  the  imposed  volume  change  AV  equal  to  the  nucleated  volumes 


from  Eq.  (33)  for  all  crack  orientation  bins.  For  the  estimate,  the 

stresses  in  each  direction  in  Eq.  (33)  must  be  approximated  in  some 

way.  Stress  appears  in  the  expression  in  two  ways:  directly,  to  describe 

the  opening  of  the  cracks  (Eq.  29),  and  in  the  nucleation  rate  function 

N*  (Eq.  31).  The  stress  appearing  in  the  nucleation  rate  function  is 

assumed  to  be  of  maximum  importance,  so  it  is  treated  more  exactly  than 

the  stress  controlling  crack  opening.  The  latter  stress,  which  appears 

directly  in  Eq.  (33),  is  taken  as  P , the  solid  pressure  at  the 

so 

previous  step.  The  pressure  is  used  here  instead  of  the  stress  in  any 
particular  direction  to  get  an  average  of  the  nucleation  behavior  in 

all  directions.  In  the  expression  for  ft*  (Eq.  31),  cr  is  treated 

C Dtp 

as  the  average  pressure  over  the  time  interval,  that  is,  (P  + P )/2, 

sn  so 

where  P is  the  nucleation-based  estimate.  The  pressure  estimate 
sn 

from  Eq.  (33)  is  then 


P 

sn 


= - P +2 a + 2a,  In 

so  no  1 


AV 

. 3 

V ft  R AtP 
CO  so 


(crack  nucleation) 
(A3) 


where  V = 


32(1 


v2) 


= 8 


A _j_\ 

^ G + C+G/3 1 


C = bulk  modulus 

G = shear  modulus 

a = nucleation  threshold 

no 

ft  = nucleation  rate  parameter 

o 

E = Young's  modulus 

v = Poisson's  ratio 

R = nucleation  size  parameter 

o 

At  = time  step 
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For  void  nucleation,  the  pressure  estimate  is  derived  from  Eqs.  (3) 

and  (4),  of  Section  II,  with  P taken  as  the  average  pressure, 

(P  + P )/2.  Then 
sn  so 


P 

sn 


= - P + 2P  + 
so  no 


2P^j0n 


AV 

3 

8-N  R At 
o o 


(void  nucleation) 


(A4) 


Usually  P is  an  overestimate  of  the  correct  pressure, 
sn 

A third  estimate  of  solid  pressure  is  made  by  considering  that 

the  imposed  change  in  volume  is  taken  by  a combination  of  expansion  and 

growth  of  existing  cracks  plus  a strain  in  the  solid  material  for  the 

brittle  case.  Equation  (37)  is  used  to  determine  the  volume  change 

associated  with  expansion  and  growth  of  cracks.  In  this  expression 

a , in  the  expansion  portion  was  replaced  by  P + cr//2,  and  a , 

CDj/j  Sg  CD0 

in  the  exponential  was  replaced  by  (P  + P )/2  + a ' where  or'  is 

sg  so 

the  most  tensile  deviator  stress.  These  replacements  for  a , were 
chosen  on  the  basis  of  trial  calculations  in  Armco  iron.  The  total 
volume  change  from  expansion  and  growth  is  then 

7 P P 

V = - T [TJ  N 1(rS3]  (P  +^)  exp  [3T  (-~  + Hr  + o'  - CT  ) At]  (A5) 
v 7“o  1 sg  2 122  go 


or  approximately 


1 3 

V = V exp  [3T  (P  + CT ' - CT  ) At]  [l  + AP  ( 7 + ^,At)]  (A6) 

v o 1 so  go  sg  p + CT_  2 1 

so  2 


where  V = the  void  volume  at  the  previous  time  step  and 
vo 


AP  = P - P 
sg  sg  so 

The  total  volume  change  is  now  taken  as  equal  to  the  combination  of 
growth  and  expansion  from  Eq.  (A6),  and  the  change  in  solid  volume  from 


74 


Eq.  17  of  this  report.  This  volume  condition  is 


AV  = V - V + V - V (A7) 

V VO  s so 

The  resulting  estimate  for  pressure  ia 


AV  - V (X  - 1)- 
vo 


V X 
VO 


1 

+ o' 


+ 


V r(E  - E )/C 
so o 

1 TE 

0 + c 

o 

2 

V 

so 

" C/po  + TE 


(crack  (A8) 
growth) 


where  X = exp  [3T  (P  + a'  - O )At]. 

1 so  go 


For  ductile  behavior  the  third  pressure  estimate  is  made  based  on  a 
combination  of  strain  in  the  solid  material  and  void  growth.  The  void 
volume  is 


V = V exp 
v vo 


P ) At 
go 


= V 


vo 


exp 


T At(P 
1 so 


T AtAP 
_1 sg 

2 


(A9) 


Combining  Eq.  17,  Eq.  (A9),  and  the  volume  condition  in  Eq.  (A7) 
provides  the  desired  estimate  of  P 

sg 
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(void  growth) 


(A10) 


AP 

sg 


AV  - V (X  - 1) 
vo 


V (E-E  ) 
so o 

c/D  + rii 

o 


V X 
vo 


T^At 

~2 


so 


C/o  + YE 
o 


where  X = exp  [t  At(P  - P )]. 

1 so  go 

The  pressure  estimate  for  the  first  iteration  is  selected  by  taking  the 
least  tensile  value  from  Eqs.  (A2),  (A3),  and  (A8)  for  the  brittle  case 
or  from  Eqs.  (A2) , (A4) , and  (A10)  for  the  ductile  case. 

After  the  first  computation  of  V and  V based  on  the  initial 

s v 

pressure  estimate,  a reestimate  is  often  required.  To  minimize  the 
opportunity  for  getting  into  loops,  two  procedures  were  constructed 
for  the  later  estimates.  The  first  is  based  on  the  result  of  the 
previous  iteration  and  on  the  physical  processes:  the  form  is  similar 

to  Eq.  (A8)  or  Eq.  (A10).  The  second  procedure  is  the  regula  falsi , an 
interpolation  based  on  two  previous  trials. 

The  first  procedure  is  based  on  a combination  of  growth,  crack 

expansion,  and  change  in  solid  volume.  After  expanding  the  exponential 

in  the  growth  function  in  the  same  manner  as  in  Eqs.  (A6)  or  (A9) , the 

crack  or  void  volume  for  a small  change  in  pressure  AP  = P - P is 

s s a 

found  to  be  related  to  the  previously  computed  volume  V , as  follows: 

va 

V = V [l  + AP  (— — T7T  + ? T,At)]  (cracks)  (All) 

v va  s P + <7/2  2 1 

a 


76 


or 


TjAt 

V = V [l  + AP  ] (voids)  (A12) 

v va  s 2 


where  P is  the  pressure  computed  in  the  previous  iteration.  The 
a 

change  in  solid  volume  is  related  approximately  to  the  change  in  pressure 
as  follows 


V 

s 


+ 


ns  \ rPs 

C ) C + P 


(A13) 


The  total  volume  change  computed  in  the  preceding  iteration  was 


AV  = V -V  + V -V  (A14) 

a va  vo  sa  so 

where  the  subscript  "a"  refers  to  the  previous  iterations,  and  "o" 
refers  to  conditions  from  the  last  time  step.  Equation  (A14)  is  sub- 
tracted from  Eq.  (A7)  to  obtain  the  change  in  AV  required  to  produce 
a correct  result. 


AV  - AV  = V - V +V-V  (A15) 

a v va  s sa 

When  Eqs.  (All)  or  (A12)  and  (A13)  are  Inserted  in  Eq,  (A15)  and  solved 

for  AP  t we  obtain 
s 


AP 

s 


V 

va 


( 


AV  - AV 

a 


1 ^ 3 T / 

Ul  * il) 

1 

P + O'/ 2 2 TlAt) 

P C 

| C + Pa 

a / \ 

0 ° 1 

(cracks) 


(A16) 
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or 


AV  - AV 

AP  = Z~T1 (voids)  (A17) 

1 

C + P 

a 

The  second  reestimate  procedure  requires  information  from  two 

previous  iterations:  the  estimated  changes  in  solid  volume,  AV  and 

sa 

AV  . and  the  computed  total  volume  change,  AV  and  AV  . The 
sb  a b 

subscript  b refers  to  values  saved  from  some  earlier  iteration.  The 
next  estimate  is  then  simply 


va  2 


AV 

s 


AV  - AV 
sb  sa 

AV  + 

sa  AV  AV 

b a 


(AV  - AV  ) 
a 


(A18) 
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Appendix  B 


INSERTION  OF  BFRACT  AND  DFRACT  INTO  HEMP 

During  brittle  fracture  calculations,  BFRACT  and  DFRACT  replaces  the 
usual  routines  for  computing  pressure  and  deviator  stress.  This  section 
describes  the  changes  required  to  incorporate  BFRACT  and  DFRACT  into 
HEMP  or  other  two-dimensional  Lagrangian  wave  propagation  codes  and 
gives  test  cases  for  verifying  the  code  results. 

Changes  to  HEMP 


The  insertion  of  BFRACT  and  DFRACT  requires  some  added  COMMON 
storage,  additional  reads  for  property  data,  a procedure  for  switching 
to  the  fracture  routine,  a CALL  statement,  and  a means  for  printing  the 
computed  damage.  These  changes  were  all  made  to  FIBROUS,  a comparable 
two-dimensional  program  at  SRI.  The  additional  COMMON  quantities  are: 


TSR  (6,30) 


ENM  (K, J) 
ENT  (K,J) 

LS 


NFR  (M) 


An  array  containing  the  fractures  and  frag- 
mentation parameters.  It  now  provides  for 
6 materials  and  30  parameters  each,  although 
only  9 parameters  are  now  in  use. 

Array  for  the  fraction  of  crack  oi  void  volume. 

Array  for  the  number  of  cracks  or  voids  per 
unit  volume  in  each  cell. 

Initializing  indicator.  Set  to  zero  at 
beginning  of  program,  reset  to  1 in  BFRACT  or 
DFRACT  after  initialization. 

An  indicator  array  for  the  type  of  fracture 
considered.  NFR  = 1 for  ductile,  and  NFR  = 2 
for  brittle  fracture. 


The  additional  materials  data  are  inserted  in  the  initializing 
routine  with  the  rest  of  the  materials  data.  In  FIBROUS,  the  indicator 
NFR  is  read  with  the  material  name.  If  NFR  equals  2,  then  two  data 
cards  are  road  to  obtain  the  first  9 variables  of  the  TSR  array.  For 
NFR  = 1,  only  one  data  card  is  needed. 
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BFRACT  and  DFRACT  are  not  called  or  Initialized  until  the  tensile 
stress  in  some  tell  exceeds  the  threshold  stress  for  nucleation,  TSR 
(Mf5).  On  the  first  CALL  (described  in  detail  below)  the  dimensional 
arrays  in  BFRACT  are  zeroed,  and  several  coefficients  that  depend  on 
material  properties  are  computed.  If  a second  material  is  involved  in 
fracturing,  these  coeffici tints  are  computed  on  the  first  call  to  BFRACT 
for  a cell  of  that  material.  DFRACT  requires  no  initializing.  As  many 
as  SuX  materials  may  undergo  f 'acture  at  once  with  the  present  array 
dimensions. 

The  CALLs  to  BFRACT  and  DFRACT  are  inserted  in  the  subroutine  that 
controls  stress  calculations  (VQP  in  HEMP,  SWEEP?  in  FIBROUS)  just  before 
the  pressure  is  computed.  The  following  statments  may  be  used: 

480  IF  (NFR(M)  - 1)  600,  480,  490 
IF  (ENT(K, J).GT.O.)  GO  TO  500 
IF  (P(K,  J)  .LT.TSR  (M.5))  GO  TO  500 
GO  TO  600 

490  IF  (ENT(K, J)  .GT.O.)  GO  TO  520 

IF  (AMIN  1 (TXX(K,  J) , TYY  (K,  J),  TTT  (K,J))  .LT.  TSR  (M,5)) 

GO  TO  520 
GO  TO  600 

500  CALL  DFRACT  (SXXN,  SYYN,  STTN,  SXYN,  EXXH,  EYYH,  ETTH,  EXYH, 
P'K.J),  ENM(K,J),  ENT(K, J) , VW,  VN,  DELTH,  E(K,J),  EE.T, 

CA(M)  , EQSTG(M),  MU(M) , RHO(M) , TSR,  YY (M)  , YD(M) , F,  M, 

DROT) 

GO  TO  540 

520  CALL  BFRACT  (LS,  SXXN,  SYYN,  STTN,  SXYN,  EXXH,  EYYH,  ETTH, 

EXYH,  P(K, J) , ENM(K, J) , ENT(K, J) , VW,  VN,  DELTH,  E(K,J), 

EEST,  CA(M) , EQSTG(M),  MU(M) , TSR,  YY (M) , YD(M) , F,  K,  J, 

M,  I CYCLE,  RHO(M),  DROT) 

540  SXXW  = SXXN 
SYYW  = SYYN 
STTW  - STTN 
SXYW  = SXYN 

600  GO  TO  650 

(usual  relations  for  pressure  and  energy) 
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The  parameters  in  the  preceding  statements  are  defined  in  the 
Nomenclature  list  given  below.  It  should  be  noted  that  BFRACT  and 
DFRACT  compute  both  deviator  stress  and  pressure.  Both  of  these 
quantities  are  positive  in  compression. 

The  damage  information  for  all  cells  undergoing  brittle  fracture 
is  listed  in  a separate  CAIL  to  BFRACT.  The  form  of  the  CALL  is: 

CALL  BFRACT  (2) 

This  statement  may  be  inserted  in  the  stress  controlling  subroutine  or 
in  a printing  or  editing  subroutine.  An  example  of  the  listing  given 
is  shown  later.  The  damage  information  computed  in  DFRACT  is  all 
contained  in  the  NM  and  NT  arrays,  which  are  parameters  in  the  CALL 
statement  and  are  available  to  the  calling  program.  Therefore,  for 
ductile  fracture  the  damage  listing  is  performed  outside  DFRACT. 

To  aid  the  user  in  testing  the  BFRACT  and  DFRACT  subroutines, 
sample  tapered-f Iyer  impact  calculations  are  given.  The  planar,  two- 
dimensional  geometry  of  the  target  and  flyer  is  shown  in  Figure  B-l. 

Sample  results  from  various  stages  in  the  two  calculations  are  shown 
in  Figures  B-2  through  B-10.  No  motion  is  permitted  in  the  third 
dimension,  and  no  vertical  motion  is  permitted  along  the  upper  and  lower 
boundaries  of  either  target  or  flyer.  The  left  and  right  boundaries 
are  free  surfaces.  The  interface  surfaces  of  target  and  flyer  are  in 
contact  during  the  early  part  of  the  calculation  but  are  allowed  to 
separate  gradually  when  the  stresses  in  the  adjacent  cells  become 
tensile.  Specific  information  about  the  calculations  is  contained  in 
the  listing  of  the  INPUT  data  decks,  Figures  B-2  and  B-7.  The  nomenclature 
for  the  data  decks  is  given  in  the  Nomenclature  list.  The  data  in 
Figure  B-2  do  not  correspond  entirely  with  the  set  of  parameters 
finally  chosen  as  the  best  representation  of  XAR  30  steel. 
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The  sample  printouts  from  the  flyer  computations  include  the 
listing  of  the  input  data,  Figures  B-2  and  B-7,  a portion  of  edits 

printed  in  SWEEPF  (Figures  B-3  and  B-8).  At  this  time  the  target  has 

passed  through  a tension  and  fracturing  phase  and  is  now  undergoing 
recompression.  Figures  B-4  and  B-9  list  a summary  of  the  damage  quantities 
at  each  cell  (printed  in  SWEEPF).  The  detailed  listing  of  the  damage 
quantities  given  in  Figure  B-5  is  printed  by  BFRACT.  The  CL  and  CN 
quantities  at  each  cell  are  given  with  the  average  crack  radius  parameter 
(CL-AV<~)  and  total  number  of  cracks  (CN-TOT).  The  total  crack  area  is 
given  and  the  cumulative  rotation  (ROT)  in  radius  is  provided.  The 
quantity  ROT  gives  the  rotation  of  the  crack  bins  from  their  original 

orientation;  it  is  not  necessarily  the  total  cell  rotation.  Figure  B-6 

contains  a portion  of  the  historical  listing  of  some  stresses.  From 
left  to  right  the  quantities  are  TXX<3,8),  P(3,8),  TXX(3,9),  P(3,9), 

TXX(3 , 10) , P(3 , 10) , TXX(4,2),  P(4,2),  P(4,2),  TXX(4,3),  and  P(4,3). 

Figure  B-10  contains  a similar  historical  listing  for  stresses  and 
pressures  at  coordinate  points  (4,4),  (4,5),  (4,6),  (4,7)  and  (4,8). 

The  computed  fracture  and  fragmentation  damage  is  shown  schematically 
in  Fig.  IV-1  of  Section  IV. 

The  fragmentation  subroutines  for  DFRACT  and  BFRACT  are  listed 
at  the  end  of  this  appendix. 
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NOMENCLATURE  OF  BFRACT,  DFRACT,  AND  THE  CALLING  PROGRAM 


Formal  and  External  Parameters 


DROT 

DTO  (or  DELTH) 

E 

EEST 

EQSTCM  (or  CA) 
EQSTGM 

ELMU  (or  MU) 

EXX,  EYY,  ETT,  EXY 

F 

J 

K 

LS 


M 

NFR 


NM  (or  ENM) 

NN  (or  I CYCLE) 
NT  (or  ENT) 

P 

RHOS  (or  RHO) 


Cell  rotation  during  time  increment, 
positive  counter-clockwise,  in  radians 

Time  increment,  sec 

Internal  energy  at  beginning  or  end 
of  time  step,  erg/g 

Estimated  internal  energy  based  on 
constant  P through  time  step,  erg/g 
, 2 

Bulk  modulus,  dyn/cm 
Gruneisen  ratio 

, 2 

Shear  modulus,  dyn/cm 

Strain  increment  in  the  x,  y,  and  8 
directions,  shear  strain  (f  + r ) 
in  the  XY  plane  xy  yx 

Thermal  strength  reduction  function 

Lagrangian  coordinate  in  the  Y direction 

Lagrangian  coordinate  in  the  X direction 

Initializing  indicator 

0 initialize  on  this  CALL 

1 computations  only 

2 print  only 

Material  number 

Indicator  array  for  type  of  fracture 

0 no  fracture  model 

1 ductile  fracture 

2 brittle  fracture  and  fragmentation 
Relative  crack  or  void  volume 

Time  increment  number 

, 3 

Crack  or  void  density,  number/cm 
, 2 

Pressure,  dyn/cm 

, 3 

Initial  solid  density,  g/cm 


SXX,  SYY,  STT,  TXY 


TSR  (1) 

TSR  (2) 

TSR  (3) 

TSR  (4)  and  (6) 


TSR  (5) 

TSR  (7),  TSR  (8) 
TSR  (9) 


TXX , TYY,  TTT 


Deviator  stress  in  the  X,Y,  and  0 
directions,  and  shear  in  the  XY  plane, 
dyn/cm^ 

(With  an  appended  N or  EN,  the  quantity 
pertains  to  the  previous  time  step; 
with  an  appended  W,  it  refers  to  the 
end  of  the  current  time  step.) 

2 

Growth  constant  = 3/(4*ETA),  cm  /dyn/sec 

2 

Growth  threshold,  dyn/cm 
Nucleation  radius  parameter,  cm 
Parameters  in  the  nucleation  function, 

o n 

no. /cm  /sec  and  dyn/cm 
N = TSR (4)*  EXP ( (P—TSR(5)  )/TSR(6)  ) 

2 

Nucleation  threshold,  dyn/cm 
Not  used 

Switch  to  indicate  whether  stress  or 
deviator  stress  governs  nucleation 

0 stress 

1 deviator  stress 

Total  stress  in  the  X,Y,  and  0 
directions 


VO  (or  VW),  VOLD  (or  VN) 

Y (or  YY) 

YD 


Relative  volume  at  end  and  beginning 
of  time  step 

, 2 

Yield  strength,  dyn/cm 

2 3 

Work  hardening  modulus,  dyn/cm  (g/cm  ) 


Internal  Variables 


CL 

3 

Cube  of  crack  radius  parameter,  cm 

CN 

3 

Crack  density,  number/cm 

DELV,  DELVA 

Imposed  and  computed  total  volume 
change  cm3/g 

DOLD,  DH 

Density  at  beginning  or  end  of  time 
step,  g/cia3 

DPJ 

Permitted  step  in  pressure  used  in 
iteration  control,  dyn/cm3 

DV,  DVO 

O 

Imposed  change  in  specific  volume,  cm  /g 
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4 


DVSA 

FNUC 


NANG 

NLOOP 

PG 

PJ 

PN 

PS 

RED.RED1 

SDH 

VSO 

WO.VV.WA 


Change  in  solid  volume,  cnrVg 

Fraction  of  cracks  nucleated  for  each 
crack  orientation  group;  FNUC  equals 
solid  angle  subtended  by  the  group, 
divided  by  4n 

Number  of  angular  orientation  groups 

Number  of  subcycling  iterations  required 

Pressure  estimate  based  on  strain, 
growth  and  expansion,  dyn/cm2 

4 

Pressure  estimate  selected  to  start 
an  iteration,  dyn/cm^ 

Pressure  estimate  based  on  nucleation 
dyn/cm2 

Pressure  estimate  based  on  strain, 
dyn/cm2 

Damage-related  reduction  factors 
for  deviator  stress 

Maximum  value  of  deviator  stress, 
dyn/cm2 

Solid  volume  at  beginning  of  time  step, 
cm3/g 

3 

Total  crack  or  void  volume,  cm  /g 
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Input  Variables 


(listed  in  order  of  occurrence  in  the  input  deck) 


I MAX 

Maximum  number  of  time  steps  permitted 

MIRR 

Indicator;  9 means  tapered  flyer 

NMTRLS 

Number  of  materials 

IJBUND 

Geometry  type;  3 means  planar  geometry 
with  laterally  immovable  edges 

I FCUT , T I MW  ITH , D ELTAT 

Unused 

COSQ 

Coefficient  of  quadrative  artificial 
viscosity  term 

XD(2) 

Flyer  velocity,  cm/sec 

TS 

Stop  time,  sec 

TRIQ 

Coefficient  of  "triangle  Q"  artificial 
viscosity  for  distortion  control 

I PRINT,  KPRINT 

Print  frequency  indicators 

NKED,  NJED 

Numbers  of  historical  edits  requested. 

Following  these  indicators  is  a large 
array  with  sets  of  4 numbers  each.  In 
each  set  the  first  number  indicates 
target  (1)  or  flyer  (2).  The  second 
indicates  the  variable  requested: 

TXX(l)  or  P(2).  The  third  and  fourth 
are  K and  J. 

3 

Original  density,  g/cm 

A three-digit  indicator.  The  tens 
column  gives  the  fracture  indicator, 

NFR. 

A three-digit  indicator  for  deviation 
stress  and  pressure  parameters.  Here 
it  indicates  that  TENS  will  be  inserted. 

Fracture  parameters 


RHO 

CFP 

DPY 

TSR 
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CA 


CB,  CC 
EQSTE 
BQSTG 
BQSTH 

PMIN 

YCC,  YCT 

TENS 

NZONES 

DEIjX 

DELY 

KROWS 

NMAT 

NCELLS 

MATR 

ITAPERED 

XLTOP 


Bulk  modulus  (dyn/cm2)  and  first 
coefficient  of  the  Hugoniot  expansion 

2 3 

P = CA  • u + CB  • u + CC  • p 

where  u = DH/RHO  - 1. 

2 

See  CA,  dyn/cm 
Sublimation  energy,  erg/g 
Gruneisen  ratio 

V - 1,  where  y is  the  exponent 
describing  expansion  of  a polytropic  gas 

Maximum  tensile  stress  permitted, 
dyn/cm2 

Yield  strengths  in  compression  and 

O 

tension,  dyn/cm 

Spall  strength  at  the  interface, 
dyn/cm2 

Number  of  zones  with  identical  K-rows 
in  either  target  or  flyer  (target  is 
treated  first). 

Cell  dimension  in  the  X direction,  cm 

Cell  dimension  in  the  Y direction,  cm 
(same  for  target  and  flyer) 

Number  of  K-rows  (or  cells  in  the  X- 
direction)  within  the  zone. 

Number  of  different  materials  in  the 
Y-direction  in  each  K-row. 

Number  of  cells  in  the  first  material 
in  the  Y-direction 

Material  number  for  the  first  cells  in 
the  Y-direction 

Indicator  for  a tapered  flyer 

Thickness  of  flyer  at  the  top,  cm 

Thickness  of  flyer  at  the  bottom,  cm 


XLBOT 


I 


target 


38. 1 mm  ■ 


(bl  ARMCO  IRON  FLYER  AND  TARGET 


MP-2024-37 


FIGURE  B-1  PLANAR  TWO-DIMENSIONAL  GRIDS  FOR  CALCULATING  TAPERED-FLYER 
IMPACTS  IN  1145  ALUMINUM  AND  ARMCO  IRON 
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FIGURE  B-4  LISTING  OF  THE  VOID  VOLUME  AND  NUMBER  OF  VOIDS  IN  EACH  CELL 
IN  THE  SAMPLE  TAPERED-FLYER  IMPACT  IN  ALUMINUM 
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FIGURE  B-5  PARTIAL  LISTING  OF  STRESS  HISTORIES  IN  THE  SAMPLE  TAPERED-FLYER 
IMPACT  IN  ALUMINUM 
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FIGURE  B-7  PARTIAL  LISTING  OF  STRESSES  AND  POSITIONS  IN  ALL  CELLS  AT  2.744 
Msec  AFTER  THE  TAPERED-FLYER  IMPACTED  THE  ARMCO  IRON  TARGET 
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FIGURE  B-8  LISTING  OF  THE  CRACK  OPENING  VOLUME  AND  NUMBER  OF  CRACKS  IN 
EACH  CELL  IN  THE  SAMPLE  TAPERED-FLYER  IMPACT  IN  ARMCO  IRON 
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1 .501)6  *06 

9,  n37E«l;S 

0* 

7 . 7 1 *#£  OS 

u. 

CN-TOT 

3.1766*06 

HOT  a 

3, JSyt-n» 

CELL 

CL  ■ 

ft.5blt-|i« 

2.2Slt-l  « 

0* 

1.873E-"* 

u. 

CL-ftvG 

3.6896-04 

p|«N*N**2  ■ 

1 . ly»L«00 

6 3 

CN  ■ 

1.39bE*(.ft 

B . 5.b£ • y 5 

0* 

6.»U*E«05 

0. 

CN-TOT 

2.6906*06 

HOT  a 

6.896E-0* 

CELL 

CL  • 

ft.b24£«04 

2.217E-V,* 

0* 

l.SS«C»P* 

0. 

CL-AVG 

3.6576-04 

PJ4N«M««2  ■ 

1 • U69E ♦ 08 

6 4 

CN  ■ 

i,3ft4E«r,6 

8.5blt«(,5 

0* 

*.»*»EA8S 

0. 

CN-TOT 

2.8596*06 

HOT  ■ 

9.95.£.o» 

CELL 

CL  ■ 

4, 36  ?t-p4 

2.191t-i« 

0* 

laSluE-O* 

u. 

CL-AVG 

3.5636-04 

P 1 •N*H**2  • 

9.73BE-IJI 

6 5 

Cn  » 

1.3376*6* 

8a23S£*l'S 

0* 

S.71oE*05 

1. 

Cn-TOT 

2.7326*06 

ROT  ■ 

1.228E-P3 

CELL 

CL  • 

ft.  1 IbE-pft 

2.10SE-0* 

0* 

la«3»E-04 

0. 

CL-AVG 

3.3ft6£-0ft 

PJ*NaR4#2  ■ 

8.S30E-01 

6 6 

CN  • 

1 .319£*r>6 

».222£«nS 

0* 

5.63IE»05 

0. 

CN-TOT 

2.7u5E*06 

ROT  a 

1.37SE-03 

CELL 

CL  ■ 

3*7ftftfc«0ft 

l,9b2£-n* 

0* 

1.3S2C-0A 

0. 

CL-AVG 

3.O806-0A 

P I •N*H**2  • 

7.039E-01 

6 7 

CN  a 

1.28bE*ob 

».13*E«0S 

0* 

S.l29£«ni 

u. 

CN-TOT 

2.6156*06 

HOT  ■ 

1 . 433E-03 

CELL 

CL  • 

3.41  QE-(«ft 

1.85.E-0* 

0. 

1.208E-0* 

L. 

CL-AVG 

2. 793E-0ft 

P1#N#H**2  ■ 

S..S1E-01 

ft  b 

CN  ■ 

1 *?u7E*oft 

7a57SE*05 

0* 

♦ a92(j£*0S 

0. 

CN-TOT 

2.4566*06 

ROT  ■ 

1.29SE-03 

CELL 

CL  ■ 

3. lulE-lft 

1.8.3E-0. 

0 « 

la!3«E-0* 

(it 

CL-AVG 

2.5636-04 

PI«N«H442  a 

A.2SSE-01 

6 9 

CN  ■ 

1 • 16  IE*  Oft 

6.T28E*<i£ 

0 • 

Aa391E»0S 

0. 

CN-TOT 

2.2736*06 

HOT  a 

»a3*2t*0» 

CELL 

CL  ■ 

2.736E-04 

1.262E-0* 

9 1 

lal*7£>«« 

0. 

CL-AV0 

2 .2656-04 

P I *N4R4#2  ■ 

2.94UE-01 

6 10 

CN  a 

l*07lE«0ft 

Sa0l8E*0S 

0 • 

*.598£.0S 

0. 

CN-TOT 

2.0346*06 

ROT  • 

2.A70E-0* 

CELL 

CL  i 

3#614E-oft 

2.S86E.0* 

*..T»E-05 

2.«»9£.04 

0. 

CL-AVO 

ft. 3956-04 

PI*N*fl»*2  ■ 

1 ,831E*00 

7 2 

CN  ■ 

1.529£*06 

«.JTSE«C5 

1 aB78E*05 

8.166E«n5 

0« 

CN-TOT 

3.471E«06 

ROT  ■ 

2.709E-0* 

CELL 

CL  ■ 

3*5l7E-oft 

S.48AE-0A 

0 • 

la93SE-0* 

o. 

CL-AVO 

♦.416E-04 

RI#N«Rft*2  • 

1..9SEA00 

7 3 

CN  • 

1 .4966*06 

*.il*(AO* 

8 • 

T.2SSE««S 

u. 

CN-TOT 

3.  lft3E*06 

ROT  ■ 

S.901E-0* 

CELL 

CL  • 

♦.8276-04 

2.UtC-0« 

0 * 

la»»7£-C4 

0. 

CL-AVG 

3.719E-04 

PI«N*R4*2  • 

l.llSlAOO 

7 ft 

CN  • 

1.387E*q6 

I.595E.05 

0* 

A.SIlOgS 

Oj 

CN-TOT 

2,8996*06 

ROT  ■ 

1.21*1-04 

CELL 

CL  • 

ft. 3I1E. oft 

t.tm-o. 

0. 

1. »!*<••* 

0. 

CL-AVG 

3.512E-06 

P!*N*Rft*2  ■ 

*.*.((•01 

7 5 

CN  ■ 

1 *3i9E«oft 

•alO*(«OS 

0 • 

•l 

CN-TOT 

2.8046*04 

ROT  • 

1.023E-03 

FIGURE  Ej-9  DETAILED  LISTING  OF  FRACTURE  DAMAGE  IN  EACH  CELL  IN  THE  SAMPLE 
TAPERED-FLYER  IMPACT  IN  ARMCO  IRON 
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ODD  OOOOOOOOOOOO  OOOO 


Subroutine  DFRACT 


SUHPQUT  If  F.  UEwaCT  (S-X,;>YY,STT.TAY.fXXI  .EYY1  ,ETT  l .f  X Y 1 .P.NM.NI  .VO.  DFR? 


I VOLO.l'fO.tiEEST  .EuSTCM,£gSTGM,ELMU,RHOS.TSR*Y,YD.F  .m.ALFAJ  OFH? 

OER? 

PJ  ESTIMATE  UE  PRESSURE  DFR? 

p«  COMPUTED  PRESSURE  BASED  ON  PJ  UFR? 

PN.PG  PRESSURES  ASSOCIATED  »ITh  NUCLEaTION  AND  GROWTH  OFR? 

MM  RELATIVE  VOID  volume  dfr? 

NT  VOID  DENSITY,  NUMHER/CM3  DFR? 

T SR ( 1 ) GROWTH  CONSTANT  * 3/(a»ETA)  UFR? 

T SR ( 2 ) GROWTH  THRESHOLD,  0YN/CM2  DFR? 

TSR  ( 3 ) NUCLEATION  RADIUS  PARAMETER,  CM  DFR? 

tsr(*i,  parameters  in  the  nucleation  function  i dfr? 

TSR(J)  i, DOT  ■ T4«EXP( (P-TSRI5)  >/TSR(6>  ) IJFR2 

T SR  ( c / NUCLEATION  ThREShoLU,  UYN/CM2  DFR? 

WO,  Vv A y/OIO  VCL'.Mt,  CM3/G  DFR? 

VGA  VOID  VOLUME  ASSOCIATED  WITH  GROWTH,  CM3/G  DFR? 

t/N a VOID  VOLUME  ASSOCIATED  wITH  NUCLEATION,  CM3?G  DERx 

DFR? 

DIMENSION  TSh ( N , 3o ) DFR? 

REAL  NM,nT,mun  DFR? 

DATA  SMF/j.Hb/  DFR? 

IF  (NM  ,LT.  i.I  RETURN  DFR? 

I OLO*RHOS/VUL(1  dfr? 

v tfo*NM'*vOLi>/Rnns  DFR? 

VSO*VOLO/RhO‘j-VVO  dfr? 

f SO*P/ (VS0*oOLfil  S PCLD*R  DFR2 

DVn*DV«(VU-VULO)/HHOS  dfr? 

IF  ( TSR (m, 7 1 .EO.  n.)  TSR(M,7)««."3.I*16*T3R(M,3I#*3#TSR(M,4)  DFR? 

•••••  DFR? 

HEOlN  SUflCVClING  LOOP  FOR  CASE  OF  URGE  STRAIN  DFR? 

«••••  DFR? 

NLOflPwMAXI  ( 1 .,-?.«DV»EuSTcM/VS0/TSRIM,5) ♦ J,5,?.5*TSR Im, 1 )*DTO*  DFR? 

I A'-IM  (P-TSR(M,2)  ,TSR(M,?»  > ) DFR? 

10P  PEI VwDW/NLOOP  DFR? 

EXX»EXX|»DELY/OVO  * ETY»ETYi»OELV/l)VO  DFR? 

ETT«ETTI»UlLV/|)VO  % £AYwEXYl*OEt.V/DVO  OFH? 

VHsVOlD/RhL'S  DFR? 

Y T ■ Y DFR2 

EO«F  DFR? 

DE«(EEST-E)?NLOOP  DFR? 

DT«nELV/OVO*DTO  DFH2 

AlsTSPCM, I ) *OT  DFR? 

DPU«P.2«(AbS(  TSR  (M,  <5  I I *ABS(P) ) DFR? 

DO  38n  NL=l«NLOGP  DFR? 

VHcVH.OELV  % El«tO*DE  DFR? 

DHwI./VH  OFH? 

TEMPI  ■ 1 ,/HMOS*fOSTGM«EI/E(JSTCM  DFR? 

c DFR? 

C ESTIMATE  OF  PRESSURE  RASED  ON  STRAIN,  GROWTH,  NUCLEATION  DFR? 

PS*PG»PN«£USTcm*( TEMPI? (VS0*0ELV)-1 .)  DFR? 

IF  (DEUV  ,GT.  <*.)  PN«?.#TSR(M,ft)»»L0G(0ELV*DH?TSR(M,7)/DT)*  DFR? 

1 2,*TSR(M.5)  - PSO  DFR? 

IE  (VV0  .IE.  0.)  GO  TO  ISP  DFR? 

XN«r.  * XP«I,0  DFR? 

IF  (PSO  .lT.  TSR ( M, S ) ) XNaTSR(M,T)/DM»oT»EXP( (PS0.TS»(M,5) )/TSR(M,DFR? 
|H1I  DFR? 

IF  (PSO  .LT.  TSR (M,?) ) XP«EXP(»i«(PS0-TSR(M»2) ) ) DFR? 

PG  * PSO* (UELV-WO*  (XP-1  . ) -XN-EOSTGM* (El-EP) /EGSTCM) ? ( VV0*XP*A1/?,DFR? 

I *XN/2,/TSH(m.6)-1 ./HhOS/EQSTCm)  DFR? 


IF  (PG  .61.  TSR(M,?))  PG  ■ PSO  ♦ (DELV-WO*(XP-| . )-F0STGM*(El-Enl/DFR2 


3 
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10 
1 1 
1? 

13 
U 
IS 
lh 
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Subroutine  DFRACT  (Continued) 


1 EOSTCm)  *E(.STCM*RHOS 

2 F (DELV  .Of.  n.  » ANU  • ►'SO  .LT.  TS«(M,2))  PG  = AM  I N ) ( PG , TSR < M, 2 ) ) 

1 Sn  PJ*ama*i (PS.PG.PN) 

Pj  * PSO  ♦ SIGNUMIMl  (AdS(Pj-PSO)  «10.«DPj>  *Pj-PSO> 

UYSsTFMP  l / ( I..PJ/EUS1CM)  -VSC 

VvA*WO*UtLV-nvS 

NC  = " . 

l;  ••  !*•••• 

C BEGIN  HEPaTIU.m  loop 



200  NL«NC«1 

VV*uVO»DELV-L)VS 

Pa  s fjSTCM* ( TfMPj/ ( vSO*DvS) - i . ) 

PNsAMJNl  (Ij.S#(PA»PSC)-TSP(M*5)  «0.) 

IF  (Pn  .LT.  C’.l  PN*EXP<PN/TSF<M.h)) 

VNA*TSP(M, 7)*Pn»UT/OH 

VOAsVVO 

P G * A M I N 1 (fi.S#(PA*PSC)-TSR(Mi2)  ,0.) 

IF  <PG  .LT.  0.)  VGAsVVO«E*P( A1»PG) 

vva*vga*vna 

OVSA«TEmP I / ( 1 ..PA/EQSTCM) -VS0 

OF.LVAaUYSA*  OVa.WO 

PJ*PA 

C TEST  fop  COMPLETION  OF  ITERATIONS 

IF  (AHS(OtLVA-OELV)/VSO  .LT.  2.E-S  .AND.  AHS  ( DV S-DVS A ) / \/SO  .LT. 

I I.E-S)  GO  TO  3-0 
IF  (NC  .GE.  In)  GO  TO  460 

IF  INC  .FO.  I .OH.  uOO t NC 1 3>  »EU.  l>  GO  TO  27* 

c interpolation  to  find  nvs 

Uvs«L'VSA*  (UVSM-UVS«)  / (OELv8-0ELvA)»  (DELV-DELv* 1 
GC  TO  2*0 

2 7 0 PJ*PA. (OELV-UEI.VA)/ (vGtt»A,/2.-TEMPi/(E0STCM»PA) »VNA/2./TS»  (M,6)  ) 

PN«PJ 

IF  <VNA*OELv-OFLV A .GT.  c.  .ANU.  VNA  .GT.  0.)  P\=?.*TSF ( M . *1 • 

1 Al.OGI  <VNA»UELV-DELVA) /VNAI  ♦ PA 
Pj«AMAXl  (PJtfl.S*  (PN«PJ)  > 

PJ*PA*SI GN ( AM  1 1 1 1 ( AbS IPJ-PA)  »OPJ) *OELVA-UELV> 

OVS*TfmP1/(1.»PJ/EUSTCM)-vS0 

IF  INC  .Eg.  II  GO  TO  290 

2 fen  IF  (AHS(OELva-dELV)  .GT.  ABS ( OELVH-DELV ) ) GO  TO  2oo 

29(,  OElvBaOEL'' A 
OvSH»nvSA 
GO  TO  200 

C FNOING  HOUTIMt 

C 

3(1  j NM«WA*UM 

NT«NT*DM/i)0L0»TSP  (M«  a)  #PN*l)T 
IF  (Nm  .GT . 0.6>  GO  TO  **0P 
bt  T *■?  . *T  x V • ALF  A/NLOOP 
ELMI'F«?.*tLMU*AMAXI  ( 1 . -SMF*VVA»UM.O.  ) 

«5l»0.666T#(OOLO“OH)/( HOLD* OH) 

UY»TXY-ELMUF/?.*E*Y« (SYY-SXX) *ALFA/NLOOP 
SaX*sx*-ELF|UF«(EXX-*<S1  ) ♦bET A 
SYYaSYY-ELMUF»(ETY-«Sl ) -BETA 
STTsSTT-ELMOF*(ETT-»(S1  ) 
aS**SXX**2*SYY**2*STT**2*?.#TXY#*2 
YE=Y»F«AMAX1 ( I ,-a.*VVA«OH,0. ) 

IF  (WS4  .LT.  YF**2/1.5)  GO  TO  3*0 
*S3aYE/’S(jRT  ( 1 .5*»S*> 

PTFOM.  (OOLO-LH)  / (OOLO*L/M)  /OT/TSR(P»  1) 


OFF?  62 
OFW?  63 
UF  H?  64 
DF»?  66 
OFM?  66 
OFF?  67 
OFF?  6H 
OFF?  69 
OFF?  7w 
OFF?  71 
DFF2  7? 
OFF?  73 
OFF?  7* 
OFF?  76 
DFR?  76 
DFR?  77 
DFR?  7 R 
DFR?  7 V 
OFF?  HO 
OFF?  81 
DFR2  82 
OFF?  83 
DFR?  HA 
DFR?  86 
DFR?  06 
OFF?  87 
OFF?  08 
DFP?  09 
DFR?  90 
DFr?  91 
DFR?  92 
OFF?  93 
DFR?  94 
DFn?  95 
DFR?  96 
DFR?  97 
OFF?  90 
DFR?  99 
DFR2  100 
DFR2  101 
OFR2  102 
DFR?  103 
DFR?  104 
DFR?  105 
DFR?  106 
DFR?  107 
DFR?  108 
OFF?  109 
DFR?  llo 
DFR?  ill 
DFR?  H2 
DFR?  113 
0FR2  114 
DFR?  115 
DFR2  116 
OFF?  117 
DFR?  HR 
DFR?  1 19 
DFF?  120 
DFR?  121 
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Subroutine  DFRACT  (Concluded) 


«S5»1  ,5/TSH  (Nt  1 1/OT 

DFR2 

122 

SXX«SXX«WS3  - EXX*wS5  ♦ PTERM 

DFR? 

123 

syy«syy«wS3  - eyy**S5  ♦ pterm 

DFR? 

12* 

STT»STT«WS3  - ETT«wS5  ♦ PTERM 

DFR? 

'25 

T XYsTX Y*WS3  - EXY4WS5/2. 

DFR? 

12b 

340 

CONTINUE 

DFR2 

127 

PSO  « PJ 

DFR? 

126 

Pj  ■ PJ*(VSO»OVS) »UH 

DFR? 

129 

P»PJ 

UFR2 

131 

Y*YT 

DFR? 

132 

VVO«VVA 

DFR? 

133 

VSO»l ,/OH-VVA 

DFR? 

134 

300 

00Ln«DH 

DFR? 

135 

£«EEST* (PuLU-PJ) *UV0/2. 

DFR? 

136 

RETURN 

DFR? 

137 

C 

DFR? 

130 

C 

EN(1  WITH  separation 

DFR2 

139 

*00 

P«Y«SXX»SyY«STT*TXY»o. 

DFR? 

1*0 

NMa.ARS(NM) 

DFR? 

1*1 

RETURN 

DFR? 

142 

C 

DFR? 

1*3 

C 

PROVISION  EOH  abort  in  case  OF  ITERATION  FAILURE 

DFR? 

1*4 

*5C 

ntpy«nTmy*i 

DFR? 

1*5 

IF  (NTRY  .GE.  s)  GO  TO  46k 

DFR? 

1*6 

OV«VO/RHOS-U/[)OLO 

DFR? 

147 

NLOOPwMAXi (3.  t-4.*2.»*NTRY»OV4EOSTCM/VSO/TSR(M*5) *0.5) 

OFR? 

1*6 

GO  TO  100 

DFR? 

1*9 

46(1 

PRINT  16PC.M,S.P.0V.DEIVA,DELVH 

DFR2 

15U 

GO  TO  300 

DFR? 

151 

1600 

FORMAT  (•  ITERATION  FAILURE  IN  DFRACT t M«*I2i*  S-*E10.3* 

DFR? 

1 62 

1 • P**E 1 0 « 3 1 • 0v»«E10.3i*  0ELV«*2E10.3) 

OFR? 

153 

End  DFRACT 

DFR? 

15* 

Subroutine  BFRACT 


SUBROUTINE  dMTACT<LS*SXXEN*SYYEN.STTEN«TXYENiEXXl,EYYi,FTTl,EXYi, 

0FR2 

2 

1 R.NM,NT.VO.VOLO*UTO*E*Ef  ST,£uSTCM.E'1STGM,ElMU*TSR*Y,  Yf)*F*KS*  JS. 

bFR? 

3 

2 m,NN,RH0S*UH0T*1PRhFR) 

BFR? 

4 

c 

BFR? 

b 

c 

NEM  --  RELATIVE  VOLUME  OF  CRACKS 

HFR  2 

6 

c 

Net  --  NUMttE*  OF  Ch«CKS/UMT  volume 

BFR2 

7 

c 

T i --  CRACK  growth  COEFFICIENT.  CM2/DYN/SEC 

BFR? 

fe 

c 

t 2 — Threshold  stress  for  growth.  oyn/cm2 

BFR? 

9 

c 

T 3 --  PARAMETER  of  NUCLEAT ICN  DISTRIBUTION,  CM 

BFR? 

10 

c 

T <•  --  NUCLEATION  rate  COEFFICIENT 

BFR? 

11 

c 

T 5 — THRESHOLD  STRFSS  FOR  NUCLEATION 

BFR? 

12 

c 

Tft  --  DENOMINATOR  OF  EXPONENTIAL  STRESS  FUNCTION 

BFR? 

1 J 

c 

T 7 --  NOT  USED 

BFR? 

14 

c 

t s --  threshold  stress  for  entering  hfract 

BFR? 

lb 

c 

t«j  --  switch  to  indicate  whether  s or  sqh  governs  nuclEation 

BFR? 

16 

c 

0 STRESS  GOVERNS 

BFR? 

17 

c 

1 DEVIATOR  stress  governs 

BFR? 

18 

c 

Cn  --  CRACN  DENSITY.  NUMHER/CM3 

BFR2 

19 

c 

CL  --  CUHE  OF  CRACK  RADIUS. CM3 

BFR? 

?D 

c 

IrRhFr (6)  --  IPRHFR(I)  *1  FOR  EXTRA  BFRACT  ITERATION  PRINTOUT 

BFR? 

21 

c 

BFR? 

22 

DIMENSION  TSR(N,30) 

BFR? 

23 

DIMENSION  CL  (11  .11.5)  .CM  11*11*5)  *C0S2Th<41  .SJN2 Th 14 1 ,CL3 (5 > . 

BFR? 

24 

1 FMICISI • HOT ( 1 1 .11 ) .STh(S) * INIT  (6) 

HFR? 

2b 

DIMENSION  VCR(N)  1VCM6I 

BFR? 

26 

DIMENSION  NTRlfll.il) 

BFR? 

27 

DIMENSION  IRRHFR (6) 

BFR? 

28 

REAL  nm.nt 

BFR? 

29 

DATA  SMF.NAHO/ | .80,5/ 

BFR? 

30 

IF  !LS  .GT.  f,l  GO  TO  20 

BFR? 

31 

c 

BFR? 

32 

c 

INITIALIZATION 

BFR? 

33 

c 

BFR2 

34 

c «•* 

INITIALIZE  GENERAL  ARRAYS  - C0S2TH,  SIN?TH.  HOT*  CN.  CL*  FNUC 

BFR? 

35 

LSw  i 

BFR2 

36 

INIT(,)«INIT(2)  wlNIT{3)«INlT(4.)«IMTCS)»INIT(6l»0 

BFR? 

37 

NANGl«NANb-l 

BFR? 

38 

FM'C(|  1*0.7071  hT/NANGI 

BFR? 

39 

F NDC (NANG) «0. 292893 

BFR? 

COS?Th(1)»1.o  S SIN2TH(I)»0. 

BFR? 

41 

DO  ip  NG«2*NANGl 

BFR2 

42 

F Ni>C  (NG)  wFNuC  ( l ) 

BFR? 

43 

TwOTH«6.2H3lH53*FLOAT (NG-l 1 /FLOAT  (NANG l) 

BFR? 

44 

COS?Th<NGI «COS ( TWqTH) 

BFR? 

45 

10 

SIN?Th(NG)*SIN(TWOTh) 

BFR? 

46 

DO  15  J*1 *11 

BFR? 

47 

DO  15  Kwl.ll 

BFR2 

48 

HOT  ( K • J)  * 0 • 

BFR? 

*9 

DO  15  NG*l*NA.gG 

BFR? 

50 

lb 

CN(K,J,NG)*CL(K,J,NGlap. 

BFR? 

51 

c 

INITIALIZE  -TSR-  COEFFICIENTS  FOR  EACH  MATERIAL 

BFR? 

52 

2b 

IF  ( INIT (Ml  .EO.  Ml  GO  TO  ?5 

BFR? 

53 

TSH(M,3l»TSR(M,3)w«3 

BFR? 

54 

VCR(M)«f).*(l./ELMU*l  ./(E0STCM*ELMU/3.)  ) 

BFR? 

55 

VCN(M1»-TSR(M,3)*TSR(M,«*I 

BFR? 

56 

INIT (M| »M 

BFR? 

57 

2b 

CONTINUE 

BFR2 

58 

IF  (LS  .E(W.  2)  GO  TO  boo 

BFR? 

59 

C 

BFR? 

60 

C 

CUMPUTAT  tONS 

BFR? 

61 
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Subroutine  BFRACT  (Continued) 


C 


C *•* 

c 

c 


loo 


c 

c • •• 


121 

c 

c 


c 


c 

130 

1 

c 


1 


1 50 


C •** 

17o 

1 


11  !NM  .LI.  0.)  n£TU«N 
VV0*NM*V0LD/8h0S 
’0*VOLO/»HCS-VVO 

ovsnvo*  (VO-VOLD)  /RHOS 

D0LD*RH0S/V0L0 

PSO  = P/(VSO*UOL(i) 

R*P0T(*S,jS)  * P0*P 

set  V4LUES  pom  MULTIPLE  LOOPS  in  case  OP  large  strain 
MULTIPLE  loops  IF  STRAIN  CORRESPONDS  to  a STRESS  GREATER 
n.33*TSR(M,S) 

NLOOPsMAX 1 ( 1 . »-4.*UV*EUSTCM/vS0/TSR(M,5> ♦ U . 5 1 

OPJ*0.2* (AbSI TSW (M*S)  ) ♦ AbS ( P SO ) I 

NTRY*1 


OEl  V«OV/NLOOP 

NTPI  ( kS « JS ) ■ 1 0(,«NLOOP* 

EXX*EXX1/NLU0P*UV/UV0 

ETT*ETTI/NLOOP«OV/OVO 

VH« | ,/OOLO  s YT*Y 

OR*nELV/OVO«UPOI  S 

Al»1.«TSR(M,  l ) *U T 


NTHY 

4 £TY«ETYl/NLOOP*DV/DVO 

5 ET.Y«EXYl/NL0OP*nv/DVO 

* JE*  (EEST-EI/NLOOP  T F1*F 
UT»U£LV/OVO*DTO 

4 TEmP1*i,/RhOS»EUSTgm»E/EOSTCM 


Than 


HEGIN  -00-  LOUP  FOR  EACH  STEP  IN  STRAIN 
DO  380  M*1 « NLOOP 

yHsvw.UELV  4 UH* 1 « / vH 

El  *E  1 *0t 
TEMPO*  TEMPI 

TEMPi*i  ./HmoS«E(jSTGH«ei/EuSTCm 
S0H*AMIN1 (SXXENtSYYENtSTTEN) 

V0Pp*r . 

Do  12n  na«1«nAnG 

VOPO«VOPO«CN (K6,  JS.NA)  *CL  < K.S • JS,NA) 

VUPO*«VCR(M)*VOPO 

*****  ESTIMATE  solid  PRESSURE  TO  BEGIN  ITERATIONS  ***** 

strain  basis  for  pressure  estimate 
PS*PG*PN*EUSTcm*(TEMPi/  « VSO*OElV)-l  .) 

IF  (P  .LT.  0.)  GO  TO  130 

,-HaCK  OPENING  BASIS  FOR  PRESSURE  ESTIMATE 
PG*PSO» ( DELV -TEMPI  ♦VSO)/  1 VOP0-TEMP  J/EUST  CM) 

IF  (PG  ,GT.  0.)  PG*PS 
GO  TO  150 

NUCLEATION  BASIS  FOR  PRESSURE  ESTIMATE 
IF  (DELV  .Gt.  n.)  PN*-PSO*  2.*TSR  <M,  b>  ,?,*T5R  ( M,  e>)  * ALOO  l ABS  ( DEL  V/ 
VCR(M)/VCN(M)/DT/PSO) ) 

GROWTH,  EXPANSION,  AND  STRAIN  BASIS  FOR  PRESSURE  ESTIMATE 
XP*FXP( Ai«AMINi  (C.*PSO*SDH-TSR (M,2) ) ) 

?G*PSO*  ( DEI V-VVO«XP»V VO- VSO* TEMPO/TEMPI* VSO>  21  VVO*XP*(  1./  (PSO*SDH) 
♦Al/2.)-VSU*vS0/EUSTCM/TEMPl ) 

Pj*AMA*l (PS,PG,PN1 
DVS*TFMP1/ ( 1 .*PJ/EUSTCM)  -vSO 
COSR«CUS (2,*R) 

SINR*S IN ( 2. •«) 

COMPUTE  STRESSES  AT  TlME(N-l)  FOR  EACH  ANGULAR  GROUP 
STH(NANG) *STTEN*PSO 
00  l 7o  NA*l , NANGl 

STh(NA)»(SXr.EN*SYYEN)/2.*PS0*(SXXEN-SYYEN»/2,*«C0S2TH(NA)«C0SR- 
SIN2Th(NA)*SINR)  *TXYEN*  <SIN2Th(Na)*COSR*COS2TH(NA)*SINR) 
SINP*SIN (2«* (R«DR) ) » COSR*COS1?.*(R*DR) 1 

NC*r 

IF  ( 1PRBFH ( 1 ) ,Nt.  n GO  TO  1220 


BFk? 

6? 

BFR? 

6 3 

BFR? 

6* 

BFR? 

6S 

BFR? 

66 

BFR? 

67 

BFR? 

68 

BFR? 

69 

BFR? 

70 

BFR? 

71 

BFW2 

7? 

BFR? 

73 

BFR? 

7a 

BFR? 

76 

BFR? 

76 

BFR? 

77 

BFR? 

ia 

BFR? 

79 

BFR? 

80 

BFR? 

81 

BFR? 

82 

BFR? 

83 

BFR? 

BA 

BFR? 

85 

BFR? 

8b 

BFR? 

87 

BFR? 

88 

IIFR? 

89 

DFR2 

90 

bch? 

91 

BFR? 

92 

BFR? 

93 

8FH? 

9a 

BFR? 

95 

BFR? 

9b 

BFR? 

97 

BFR? 

98 

BFR? 

99 

BFH2 

100 

BF»2 

101 

BFR? 

102 

BFR2 

103 

BFR2 

104 

BFR2 

105 

BFR2 

106 

BFR? 

107 

BFR? 

108 

BFR? 

179 

BFR? 

111) 

BFR? 

111 

BFR? 

112 

BFR? 

113 

BFR? 

114 

BFR? 

115 

BFR? 

116 

BFR? 

117 

BFR2 

118 

BFR? 

119 

HFR? 

I2l) 

BFR? 

121 

102 


Subroutine  BFRACT  (Continued) 


HrflMT  20lO.NN,«S»JS 

RFR? 

i?? 

Ph  J \T  3101«WO 

BFR2 

123 

P><  I NT  1201  .PJ.PS.PG,PN.P«0VS.  DEL  V 

BFR? 

12o 

2l)l>  l 

CONTINUE 

BFR2 

125 

3001 

FOHMAKSXf*  WO»*E12.3) 

BFR2 

126 

1201 

FuRnaTI*  Pj=«Fll.3**«  PS**£ 1 0 . 3 • • • PG*#E 1 0 • 3 1 • • PN**Elu.3t*( 

P»*  BFR? 

127 

1 Eln.3.*.  Uvs«#£ln,3«*t  0ELv*»E10.3) 

PFR? 

I2B 

201  C 

FORMAT  (5***  NN*»K,*  K5,JS»*?2*1 

BFR? 

129 

1 22C 

continue 

BFR2 

13j 

<:*•** 

(••••BFR? 

131 

C 

BEGIN  ITERATION  LOOP 

BFR? 

132 

1 

•••••  •••< 

133 

200 

CONTINUE 

BFH? 

13* 

nC*mC» 1 

BFR? 

135 

vv*wr*i)ELV-ovs 

BFR? 

136 

c 

compute  PRESSURE 

BFR? 

137 

I 

PAaFQSTCM* ( TEMPI /( VS0*UVS )- 1.  ) 

BFR? 

1 3P 

VV=VM-VSO«DVS 

BFR? 

139 

J 

c ••• 

COMPUTE  UEVIATOR  STRESS 

BFR? 

1*0 

RED*  AmaX  l (O.i  1 .-#.*VV*DH) 

BFR? 

1*1 

RED  1 *AMAX I ( 1 .-SMF*VV*OHtO,  I 

BFR? 

1*2 

*Sl «. 66667# (UOLO-OH)/  (UOU)«OH) 

BFR? 

1*3 

btTA*?.«TXYEN*OP0T*0ELV/Uw 

BFR? 

1** 

ELMUFsREt)l*2.«ELMU 

BFR2 

1*5 

k 

Tktf*txtEiv-EI.muF/2.«ExT«  (STYEN-SXXEN»«0H0T»DELV/0V 

BFR? 

1*6 

SXXe*SXXEN-ELmuF«(EXX-»S1 1 *8ETa 

BFH? 

1*7 

i 

SYrF*SYYEN-ELMUF*(ETY-W51 1-bETA 

BFR2 

1*6 

STTE*STTEN-tLMUF*lETT-«lSl  ) 

BFR? 

1*9 

*S4*SxX£**2*SYVE#«2«STTE##2*2»*TXY£**2 

BFR? 

1 5 0 

1 

, 

YE*YT*F«REO 

BFR? 

151 

IF  (WS4  .LE.  YE«*2/1.5  1 GO  TO  230 

BFR? 

15? 

*S3*yF/50HT ( 1 .5**5* ) 

BFR? 

153 

SXXF.SXXE*1>S3 

BFR? 

15* 

SYYE*SYYE»«S3 

BFR? 

155 

TXYF*TXYE**S3 

BFH2 

156 

5TTe*STTE*»S3 

BER2 

157 

1 

23' 

CONTINUE 

bFR? 

150 

C *•• 

COMPUTATION  of  CRaCN  VOLUME  FROM  ELASTIC  opening,  gro«tm, 

BFR? 

159 

c 

NUCLEaTION  ANU  FRAGMENTATION 

BFR? 

160 

VVA»0  • 

BFR? 

161 

00  ?5n  NA*ltNANG 

BFH? 

162 

IF  (NA  .LT.  NANG)  GO  TO  237 

BFR? 

163 

1 

STHW*STTE«PA  S GO  TO  2*« 

BFR? 

164 

J 

237 

STh*«pa» (SXXt*SYYE)/2.* (SXXE-SYYE)/2.«IC0S2TH(NA)«C0SR-SIN2TM(NA)*BFR2 

165 

1 SInR)*TxyE* (SIn?TH(nA1  •COSR*CoS2TH (NA) *SlNR) 

BFR? 

166 

r 

2*0 

SAVG»(STH(NA) ♦ STh* ) / 2 . 

BFR2 

167 

i 

UTC»CN(KSt  JS»NAl*DH/UOLl>*CL<KSt  JS»NA> 

BFR? 

160 

IF  (SAVG  .IT.  TSH(M.?)I  OTC»DTC#EXP ( A 1 • ( SAVG-TSR (M*2 1 1 1 

BFR? 

169 

SCN«SAVG-TSW(M.91»(PS0*PA)/2.-TS«(M»5) 

BFR? 

170 

OTN*0. 

BFR? 

171 

‘ * 

IF  (SON  .LT.  o.)  DTN.TSR(M,A|«EXP(SCN/TSR(M,6))#DT*FNUC(NA) 

BFR? 

172 

1 *TSR(M.3I 

BFR? 

173 

1 

IF  (STH*  ,LI.  p.l  VVA»VVA-VCR<M1 "STH**! DT C*DTN> 

BFR2 

17* 

25o 

Continue 

0FR2 

175 

VVAsVVA/UH 

BFR? 

176 

c •** 

COMPUTE  CHANGES  IN  V AND  IN  V SUB  S 

BFR? 

177 

Sun* am  IN | (SXX£  tSYYEtSTTEl 

BFR2 

170 

OVSA«OVS 

BFR? 

'79 

0ELVA*0VS*VVA-VV0 

BFR? 

100 

PJ*PA 

BFR? 

181 

103 


f 


Subroutine  BFRACT  (Continued) 


IF  (1PH8FSO)  .NE,  1)  00  TO  1320  BFR2 

PRINT  1301.NC.DVS,0VSA,OElV.D£LVA,0V.VV*.PA,S0H  BFR? 

PRINT  3002.5XXE .SYYE.STTF.TXYE  BFR? 

300?  F0HMAT(5X.*  SXXEtS¥YEiSrTEtTX¥E»*<*E12.3)  BFR? 

13oi  FORMATdX,*  NC»*12»*  0VS=*2E11.3t«  DELV»*2E1 1 .3**  DV**E11.3»  BFR? 

i # VVA«*E11.3,*  PA«*E11.3f«  SDH««EU.3)  BFR? 

1320  CONTINUE  BFR2 

C BFR? 

C ***  TEST  FOR  COMPLETION  OF  ITERATIONS  BFR? 

IF  (A8S(0F.LVA-0ELV)/VS0  .LT.  2.E-51  GO  TO  300  BFR? 

IF  (NC  .GE.  in)  GO  TO  450  BFR? 

IF  (NC  .E<J.  1 .OH.  M00(Nc*1i3)  .£Q.  0)  GO  TO  270  BFR? 

C INTERPOLATION  TO  FINO  UVS  BFR2 

0VS«0VSA*  (UvSH-DVSAl/UlELVD-DELVAlMDELV-DELVA)  BFR2 

GO  TO  280  BFR? 

2To  PJ«PA* (OELV-OELVA)/(VVA«( 1 ./ (PA ♦ SDH/?.  1 ♦ A 1/2. 1 -TEMP  1? ( EQSTCM*PA 1 ) BFR? 

IF  (Pj  .LT.  0.  .OR.  PA  .GE.  0.)  GO  TO  27V  BFR? 

PJ«PA*EoSTCm«(vVA.OELV)/VSO  BFR? 

IF  (PJ  .LT.  (j.)  PJ*AMAX1  (PJ.PA) /2.  BFR? 

27V  PJ»PA*SIGN(AMINI  (AMS(PJ-PA) ,OPJ)  .UELVA-DELV)  BFR? 

DVS«TEMP1/ ( l.«PJ/EUSTCM)-VSO  BFR? 

IF  (NC  .Eg.  1)  GO  TO  290  BFR? 

280  IF  (AHS(OELYA-OELV)  .GT.  ABS (DEL VB-DELV ))  GO  TO  200  BFR? 

29o  DElvB«DELVA  * OVSBaUVSA  % GO  TO  2i,0  BFR? 

C BFR? 

C ENDING  ROUTINE  BFP2 

300  CONTINUE  BFR2 

M«r.  BFR? 

RaR.OR  BFR? 

DO  320  NA«1 .NANG  BFR? 

IF  (NA  .LT.  NANG)  GO  TO  3(j7  BFR? 

SThw«$TTE*PJ  J GO  TO  3)0  BFR? 

3C7  SThR»PJ*(SXXE*S¥¥£)/2.*(SXXE-S¥YE)/2.*(C0S?Tm(Na)*C0SR-SIN2TH(NA)*BFR? 

1 SlMR)*TX¥E«(SlN2TH(NAiACOSR*COS2TH(NA)»SlNR)  BFR2 

310  SAVG«(STH(NA).SlMB)/2.  BFR2 

5CN»SaVG-TSH(M.9) •(PS0*PJ) /2.-TSH(M*5>  BFR? 

ON*n,  BFR2 

IF  (SCN  .LT.  0.)  DN«TSH(M,*)«EXP(SCN/TSR(M,6))«DT»FnUC(NA)  BFR? 

CNO^CN (KS.JS.NA)  BFR2 

CN(kS. JS.NA)  »CN(  KS.JS.NA)  •OH/OOLO*l)N  BFR2 

IF  ( CN (KS.JS.NA ) ,EQ.  0.)  SO  TO  320  BFR? 

CL  (KS.  JS.NA)  »<CN0*CL  (KS.JS.NA)  •EXPUl*  AM  InI (SAVG-TSRlM.2) .0.) )♦  BFR? 
1 DN*TSR(M.3) )/CN(KS. JS.NA)  BFR? 

NT»nT*CN(KS»JS.MA)  BFR? 

320  CONTINUE  BFR2 

NM« VV A*DH  BFR2 

PS0»PJ  BFR2 

IF  (NM  .GT.  0.6)  GO  TO  400  BFR2 

c •••••  EVALUATE  average  GROSS  PRESSURE  FROM  SOLIO  PRESSURE  4**BFR? 
PJ»PJ*(VS04UVS)*UH  BFR? 

SXXEN»SXXE  BFR2 

S¥¥EN«S¥¥E  BFR? 

STTEN.STTE  BFR? 

rx¥FN«TX¥E  BFR? 

P«PJ  BFR? 

Y»¥T  BFH2 

TXXE«PJ*SXXE  BFR? 

IF  ( IPRBFR ( 1 ) ,Ew.  1)  PRINT  1255.NN.NL. P50.PJ.TXXE.NM  BFR? 


1255  FORMATOOM  CONVERGED. 5X.*NN»*l4.«»NL»*IA.*.  PS0»*E1 0.3.«»  PJ»*E10.BFR2 


182 

183 

184 

185 
1 86 
187 
1 88 
189 
1 9(j 
191 
) 92 

193 

194 

195 

196 
I 97 
198 
I 9V 
200 
201 
?02 
?03 

204 

?05 

206 

207 

208 
?09 
?10 
?11 
?12 
?13 

214 

215 

216 
217 
?1B 

219 

220 
221 
?22 
?23 

224 

225 

226 
?27 
228 
?29 

230 

231 
?33 

234 

235 

236 
?37 
238 
’39 
2*0 
241 


104 


Subroutine  BFRACT  (Continued) 


13.*.  TXXe*£  10.3.*.  Nh»*Elo.3) 

BFR2 

2*2 

c 

8FP? 

2*3 

c 

FNO  OF  SUBCYCLING  LOOP 

BFR2 

?♦* 

c 

BFR2 

2*5 

VVO«VVA  S V so«vso*uvs 

BFR? 

?*6 

3a  n 

00Ln«DH 

BFR? 

2*7 

E=FFST* (P0-PJ)*l)V0/2. 

BFR? 

2*8 

HOT (KS. JS) *R 

BFR2 

2*9 

RETURN 

BFP? 

250 

c 

BFR2 

251 

c 

FNO  WITH  SEPARATION 

BFR2 

252 

*00 

SXXF*SYU*S*Y£*TAYE*PwO. 

BFP? 

253 

Y*YT  * NM*-AHS(NM) 

BFR? 

25* 

return 

BFR? 

255 

c 

BFR? 

?56 

c *•« 

PROVISION  FOR  ABORT  IN  CASE  OF  ITERATION  FAILURE 

BFR2 

?57 

*50 

NTPy«nTRY»1 

BFR? 

25B 

IF  (NTRY  .GE.  5)  GO  TO  *bn 

BFR2 

259 

DVsVO/RHOS-l  ./OOLO 

BFR? 

260 

NLOOPwMAX 1 ( J..-4. *2.**NTHY*0V*£QSTCM/VS0/ T3P(M. 51 *0.5.2. *Ni.OOP> 

BFR? 

261 

NTPI (kS. jS) ■nTRY*100*NLUOP 

BFR2 

?62 

GO  TO  100 

BFh? 

263 

*6c 

PRINT  1600.NN.kS.JS. SOm .P , DV , DEL V A , DEL VB . DEL V , DvO . VO 

BFR? 

26* 

IF  (NTRY  .£(j.  s)  STOP  22 

BFR? 

265 

NTw'i, 

BFR? 

266 

H*R*OP 

BFR2 

267 

00  620  NA»l «NANO 

BFR? 

268 

IF  (Nfl  ,LT.  NANG)  GO  TO  607 

BFR? 

269 

STH,«eTTE*PJ  * GO  TO  610 

BFR2 

270 

617 

STHwwPJ. (SXXE*SYY£)/2.*(SXXE-SYYE)/2.*(COS2TH(NA)#COSR-SIN2TH(NA)*BFR2 

271 

1 CInR)*TXYE*(SIN2TH(NA)*COSH*COS2TH(NA)«SIN«> 

BFR? 

?72 

610 

5AVG»(STH(Na) «SThW)/2. 

BFR? 

?73 

ScN»SaVG-TSR(M,9)*(PS0*PJ)/2.-TSR(M.5) 

BFR? 

27* 

0 N * (■  . 

BFR? 

275 

IF  (SCN  .LT.  0.1  DNwTSR1m.*)*EXP(SCn/TSR (m.6) >*OT*FnUC(N*) 

BFR2 

276 

CNO«CN (KS.JS.NA) 

BFR? 

277 

CNIKS.  JS.NA) *CN(KSi JS.NA) *0H/00L0*UN 

BFR? 

278 

IF  (CN (KS.JS.NA)  ,EQ.  0.)  GO  TO  620 

BFR2 

279 

CL (KS.JS.NA) « (CNO*CL (KS.JS.NA) *EXP< A] • AM  INI (SA VG-TSR ( H.21 .0.) ) ♦ 

BFR? 

280 

1 0N*TSH(M,3) )/CN(KS. JS.NA) 

BFR2 

281 

NTsNT*CN(KS» JS.NA) 

BFR? 

282 

620 

CONTINUE 

BFP2 

?83 

NM.wa*OH 

BFH? 

28* 

IF  (NH  .GT.  0.6)  GO  TO  *Po 

BFR? 

285 

pjrR  i*  (uen.nvc)  «riF 

BFR? 

?86 

SXXFNaSXXt 

BFR2 

28B 

syyfn«syye 

BFP? 

289 

STTEN-STTE 

BFR? 

290 

TXYFNwTXYE 

BFR? 

291 

PaPJ 

BFR2 

292 

Y»YT 

BFR? 

?9j 

VVOwVVA  S VSOaVSO*L)VS 

BFR? 

?9* 

DOLD»DH 

BFR? 

295 

E»EEST« (PO-PJ) *DVO/2. 

BFR2 

296 

ROT(KS,JS)»R 

BFR? 

297 

return 

BFR? 

296 

c 

BFR? 

?99 

C 

final  printout 

BFR? 

300 

C 

BFR2 

301 

105 


Subroutine  EFRACT  (Concluded) 


bOO 


b JO 


blr, 


b20 
S2 1 

lbOO 

15li 

1 b 1 1 
1600 


izero=i 

PH  1 NT  1511 1 ( (NTH I {K,Ji t J»1 *1 1 I .**1.11  ) 
JO  53«  K>lill 
00  53"  J«lfll 


RFH2  102 
BFR2  103 
BFH?  30* 
BFH2  305 


MHIlK.Jln  BFH?  306 

00  521  Ksl.ll  BFH2  307 

00  52n  J*l.ll  BFH?  3o« 

IF  (CNIK.Jtll  ,Eu,  GO  TO  520  BFH?  30H 

IF  (IZEHO  .EU.  |)  PAINT  1500  BFH?  310 

1 it H0*2  BFH?  Ill 

CNSUH«CflIT2«CHlT3»CHlT**r,.  BFH?  112 

U 0 510  NA*1 • NANG  BFH?  Jl3 

CL3(NA)*CL(n.J.N»)**(1./3.)  BFH?  11* 

CHIT2«CHIT2«CN(K.JjN4)*CL3(NA)**2  BFH?  315 

CNSUM«CNSUH*CN(M J.NA)  BFH?  316 

CHIT3cC«!T3*CN<K*JiN4)*CL(K»JtNA)  BFH2  317 

CHIT2«3.1*16*CHIT2  BFH?  31H 

HAn«(CHlT3/CNSUMl»*(l./3.)  BFH?  319 

PhInT  lSlut <CL3(I> .1*1.5  .RAD.C»IT2.K.J.1CN(k.J.I).I«1.5).CNSUM.  BFH2  32o 

1 hot(K,j)  BFH?  321 

CONTINUE  BFH?  122 

CONTINUE  BFH?  323 

HETI/HN  BFH?  32* 

FOHMAT  < IHiii*  LISTING  OF  CHACO  LENGTH  AND  NUMBFH  FOH  EVE.’T  ••  FjFH?  325 

1 *FHACTUHING  CELL*/)  BFH2  326 

FORMAT  (•  CELL  CL  * **E 1 0 . 3 .2*.' 1 0 . 3 1 • CL-AVG  » *El0.3*  BFH?  327 

1 1 * h HI*N*H**2  a «E10. 3/213*  CN  ■ **E  lt<  .3.2X  *E  10 . 3,*  CN-TOT  ■ BFH?  328 

2*E 1 1 . 3i 1*H  WOT  ■ iEIO.3/1  BFp?  329 

FOHMATI*  *«10X**  HFHACT  ITERATIONS— NTRYOoo.NLI  ^*./.10X»*  K-BFR2  330 

lHO*S  ANO  J-COLS*»//.ll(llIl0/)*/*#l*»  BFR?  131 

FORMAT  <*  ITERATION  , FlLUHE  IN  BFHACT.  N**I5.*»  Ka*I3.*.  J**I3.  BFR2  332 


1 5Ei2.3,/.5x,3Ei2.J) 

End  hfract 


BFh?  333 
BFH?  33* 


106 


Appendix  C 


CALCULATION  OF  VOID  GROWTH  UNDER  HIGH  SHEAR 

The  calculations  of  this  appendix  were  made  to  guide  in  constructing 
stress-strain  relations  for  material  containing  voids  for  cases  where 
the  material  is  subjected  to  high  shear  strains.  The  material  was 
idealized  as  an  elastic-viscous-plastic  material,  initially  containing 
a uniform  spacing  of  spherical  voids.  The  "typical"  element  used  in  the 
computations  was  a cyxinder  with  height  equal  to  the  diameter  and  con- 
taining one  spherical  void  at  its  center  as  shown  in  Figure  C-l.  The 
calculations  simulated  a high  speed  extrusion  that  could  occur  in  10  to 
100  nsec.  The  planar  boundaries  were  moved  outward  to  simulate  tensile 
straining  on  the  material.  The  cylindrical  boundaries  were  moved  inward 
at  half  the  velocity  of  the  planar  boundaries  to  produce  a constant 
volume,  high-shear-strain  loading  on  the  material. 

The  computations  were  conducted  with  a computer  program  termed 
VOID,  a special  purpose,  finite-difference,  two-dimensional  wave  propaga- 
tion code  described  in  Reference  3.  In  the  calculations,  the  motion  of 
the  boundaries  caused  waves  to  propagate  through  the  body  of  the  material. 
By  use  of  a high  coefficient  of  artificial  viscosity,  the  wave  amplitudes 
were  minimized  so  they  did  not  interfere  with  the  interpretation  of  the 
steady-state  results. 

The  quantities  of  interest  for  determining  stress-strain  relations 
for  fracturing  material  are  the  pressure,  deviator  stresses,  gross 
specific  volume,  void  volume,  and  gross  shear  strain.  The  stress 
quantities  were  computed  as  averages  based  on  the  total  forces  on  the 
planar  and  cylindrical  surfaces.  The  growth  rate  of  the  void  was  derived 
from  the  computed  history  of  void  volume.  Because  only  a shear  strain 
was  imposed,  the  pressure  remained  small,  making  interpretation  of 
pressure  volume-relations  uncertain.  However,  the  void  growth  rate  does 


appear  to  follow  the  rule 


P - P 

R = — 2 R (C-l ) 

477 


where  R 
P 
P 

o 


average  void  radius 
pressure 

a threshold  pressure 


Tj  = the  material  viscosity  used  in  the  stress  relaxation 
model  for  the  material. 

Hence  the  growth  rate  is  the  same  under  high-shear-strain  loading  as  for 
spherical  or  one-dimensional  loading. 


The  deviator  stress-strain  relation  was  also  studied.  As  in  the 

3 

one-dimensional  calculations,  some  of  the  results  fit  the  form  we 
expected,  which  is,  in  one  dimension 


. , 4 V (T  - 2/3Y 

a = - G( 1 - 4V  0)  

3 v V T 


where  o' 
G 
V 

v 

P 

V 

Y 
T 


deviator  stress 

shear  modulus  for  the  solid  material 
specific  volume  of  the  void 

density 

gross  specific  volume 
yield  strength 
a time  constant. 


(C-2) 


Some  of  the  results  showed  a long  rise  time  corresponding  to  a much 
smaller  modulus  than  G.  These  latter  results  may  have  been  caused  by 
the  effect  of  waves  on  the  results.  Because  we  could  not  interpret 
these  results,  we  chose  to  use  Eq.  (C-2)  for  deviator  stress  until  a 
more  complete  study  could  be  conducted  of  deviator  stress. 
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Appendix  D 


BALLISTIC  EXPERIMENTS  ON  1145  ALUMINUM  AND  ARMCO  IRON 

As  an  aid  to  the  code  development  effort,  long  rod  ballistic  impact 
experiments  were  performed  on  two  materials  whose  dynamic  fracture 
behavior  under  one-dimensional  strain  is  well-known.  We  performed  4 
experiments  on  10  x 10  cm  plates  of  1. 27-cm-thick  1145  aluminum  and 
14  experiments  on  similar  specimens  of  Armco  iron,  using  the  testing 
facility  shown  in  Figure  14  of  Section  III.C.  The  dynamic  fracture 
parameters  for  both  materials  under  uniaxial  strain  had  been  determined 
in  previous  work.  Our  goal  was  to  count  and  measure  the  fracture  damage 
in  these  specimens  and  chen  eventually  compare  the  results  with  the 
predictions  of  the  evolving  two-dimensional  fracture  model. 

Table  D-I  summarizes  the  results.  For  the  experiments  on  armor 
steel  described  in  Section  III.C,  the  projectiles  were  right  cylinders 
1.03  cm  in  dlamter  by  2.03  cm  long,  made  of  drill  rod  and  fitted  into 
Lexan  polycarbonate  sabots.  For  the  1145  aluminum  experiments,  the 
projectiles  were  in  the  as-received  spherodized  condition  at  a hardness 
of  Rockwell  B95 . For  the  Armco  iron  experiments,  the  projectiles  were 
heat  treated  to  Rockwell  C35,  The  projectiles  were  accelerated  to 
different  velocities  to  produce  a range  in  the  extent  of  fracture  damage 
in  the  targets. 

Work  on  1145  aluminum  was  discontinued  after  four  experiments,  since 
it  became  apparent  that  no  voids  or  cracks  could  be  produced  under  these 
conditions.  The  penetration  mechanism  for  this  soft  aluminum  alloy  was 
one  of  pure  plastic  plugging,  although  spherical  ductile  voids  were 
previously  found  to  develop  under  uniaxial  strain  conditions.  The 
photomicrograph  in  Figure  D-l  of  a polished  cross  section  taken  through 
the  partially  penetrated  Specimen  No.  3 shows  the  intense  plastic  shear 
deformation. 
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FIGURE  D-1  POLISHED  CROSS  SECTION  OF  1145  ALUMINUM  PARTIALLY 
PENETRATED  BY  A ROD-LIKE  PROJECTILE 
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Fracture  damage  did  accompany  projectile  penetration  in  Armco 
iron  at  velocities  of  about  li.00  ft/sec  and  above,  as  could  be  readily 
observed  on  the  rear  surfaces,  Figure  D-2.  Full  penetration  did  not 
occur  at  velocities  up  to  2130  ft/sec. 

Examination  of  polished  sections  with  an  optical  microscope  showed 
that,  as  in  the  uniaxial  strain  situation,  the  damage  was  in  the  form 
of  planar  microcracks,  as  shown  in  Figure  D-3.  The  cracks  lay  nominally 
in  a band  about  2 mm  from  the  rear  surface  and  roughly  parallel  to  it. 
They  coalesced  most  strongly  in  the  material  directly  in  the  path  of 
the  projectile,  and  at  sufficiently  high  velocities  resulted  in  rear 
surface  spallation. 

To  obtain  the  necessary  experimental  information  for  eventual 
comparison  with  the  predictions  of  the  two-dimensional  dynamic  fracture 
model,  we  assessed  quantitatively  the  fracture  damage  in  the  specimens. 
The  procedure  entailed  counting  and  measuring  the  microcrack  traces  on 
enlarged  photomicrographs  of  polished  sections,  as  showr  in  Figure  D-4. 

A statistical  transformation  may  be  applied  to  these  aata  to  obtain 
the  actual  crack  size  distributions  per  unit  volume.  The  crack  trace 
counting  and  measuring  operation  is  facilitated  by  a large  area  record 
reader,  and  the  statistical  transformation  is  performed  by  a computer 
code.  The  counting  and  measuring  was  carried  out  for  the  above  Armco 
iron  specimens,  but  the  statistical  transformation  to  volume  distribution 
was  not  accomplished.  In  any  future  work  it  is  recommended  that  this 
be  done  to  provide  data  for  comparison  with  HEMP  code  predictions  based 
on  the  fracture  model  developed  in  this  program. 
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FIGURE  0-2  BACK  SURFACE  FRACTURE  IN  ARMCO  IRON  SPECIMENS 
AT  VARIOUS  PROJECTILE  VELOCITIES 
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1300  ft /sec 


FIGURE  D-3  POLISHED  AND  ETCHED  CROSS  SECTIONS  OF  ARMCO  IRON 
SPECIMENS  SHOWING  FRACTURE  DAMAGE  INDUCED  BY 
PROJECTILE  IMPACT 
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MP  2024-5 


FIGURE  D-4  HIGH  MAGNIFICATION  VIEW  OF  MICROCRACK  TRACES 

ON  AN  UNETCHED  POLISHED  SECTION  THROUGH  ARMCO 
IRON  SPECIMEN  NUMBER  15 


Appendix  E 


IMPROVED  HOMOGENEOUS  STEEL  ARMOR 

An  improved  homogeneous  steel  (IHS)  armor  can  be  made  by  uni- 

16 

directionally  solidifying  and  homogenizing  steel  castings.  This 
procedure  results  in  a cleaner  steel,  free  from  large  inclusions  and 
having  a lower  content  of  detrimental  impurity  elements  such  as  sulfur, 
oxygen,  and  phosphorus.  In  addition,  the  procedure  provides  a preferred 
texture  and  a fine  homogeneous  dispersion  of  included  particles.  These 
microstructural  features  are  thought  responsible  for  the  excellent  high 
strength  ductility  and  superior  ballistic  impact  properties  of  the 
material . 

The  laboratory  method  for  making  IHS  consists  of  pouring  the  steel 
under  vacuum  into  a mold  cavity  on  a thick  copper  chill  block.  The 
feasibility  of  large-scale  production  by  the  commercial  electroslag 
remelt  (ESR)  technique  is  being  actively  investigated  by  AMMRC,  and 
the  results  are  encouraging. 

Because  of  the  attractive  properties  of  the  improved  armor  steel 
and  the  promise  of  the  ESR  process  as  a cost-effective  manufacturing 
method,  some  preliminary  investigations  were  carried  out  on  IHS  in  this 
program.  The  microstructure  was  examined,  quasistatic  tensile  properties 
in  the  short  transverse  (through-the-plate-thickness)  direction  were 
determined,  the  fracture  surfaces  were  examined,  and  the  region  near 
the  hole  produced  by  a ballistic  round  was  examined. 

Material 

17 

The  steel  was  manufactured  by  U.S.  Steel  Corporation  by  the 
unidirectional  solidification  technique  in  18  x 13$  x 18  inch  molds  on 
a thick  copper  chill  block.  The  steel  had  the  following  approximate 
chemical  composition  (wt%) : 0.52  C,  0.70  Mn,  0.25  Si,  1.20  Ni,  0.75  Cr, 

0.50  Mo,  0.025  Al,  < 0.003  P,  0.005  S,  and  < 0.009  N. 
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Melting,  alloying,  and  pouring  were  done  under  vacuum.  The  vacuum 

melting  was  done  in  a 500-pound  induction  furnace;  the  pouring  tempera- 

o 

ture  was  maintained  at  about  2900  F.  Homogenization  took  place  in  a 

commercial  high-temperature  heat  treating  furnace  at  a pressure  of 

o 

1 micron  and  a temperature  of  2400  F for  64  hours.  After  furnace  cooling 
o 

to  1000  F and  then  air  cooling,  the  casting  was  cut  into  four  slabs  and 

hot  cross-rolled  to  half-inch  plates.  The  plates  were  then  heat  treated 

o 

to  about  RC  60  by  heating  to  1500  F for  1/2  hours,  oil  quenching,  and 

o 

immediately  tempering  at  250  F for  1 hour. 

Microstructural  Observations 

Polished  and  etched  surfaces  in  three  mutually  perpendicular  planes 
of  a plate  of  undirectionally  solidified  and  homogenized  steel  were 
examined  with  the  optical  microscope.  Most  notable  were  the  fineness 
and  homogeneity  of  the  microstructure,  the  absence  of  large  inclusions, 
and  the  lack  of  banding.  Figure  E-l  shows  micrographs  of  the  IHS  armor 
and,  for  comparison,  the  XAR30  armor.  Note  the  coarser  structure  and 
the  large  inclusions  in  the  XAR30  material. 

Tensile  Property  Measurements 

Tensile  properties  in  the  short  transverse  direction  (through  the 

plate  thickness)  were  determined  by  spark  machining  small  cylindrical 

blanks  about  6.35  mm  in  diameter  and  12,7  mm  long  from  the  armor  plate, 

grinding  the  central  sections  of  the  blanks  to  about  3.18  mm  to  obtain 

a typical  dumbbell  shaped  specimen  with  a gage  length  of  about  5.08  mm, 

and  then  pulling  the  specimens  in  a tensile  testing  machine  at  a constant 

crosshead  speed  of  0,5  mm/minute.  Triplicate  specimens  yielded  extremely 

reproducible  load  displacment  records,  shown  in  Figure  E-2  which  were 

analyzed  to  obtain  the  short  transverse  tensile  properties  given  in 
* 

Table  E-I.  The  curves  from  similar  tests  for  conventional  rolled 

* The  results  for  similar  specimens  tested  at  a 100-fold  faster  testing 
speed  (50  mra/minute)  were  essentially  identical  and  show  the  1 ■*  strain 
rate  sensitivity  of  the  steels  in  this  range. 
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(a)  IMPROVED  HOMOGENEOUS  STEEL 


(b)  XAR30  ARMOR  STEEL 

MP-2024-29 


FIGURE  E-1  COMPARISON  OF  CROSS  SECTIONS 
OF  IMPROVED  HOMOGENEOUS 
STEEL  ARMOR  AND  CT  XAH30 
ARMOR 
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STRESS  — 1000  Ib/ln 


RESULTS  OF  QUASI STATIC  TENSILE  TESTS 
OF  THREE  ARMOR  STEELS  TESTED  NORMAL  TO  THE  ROLLING  PLANE 
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homogeneous  armor  steels  are  included  for  comparison.  These  properties 

are  of  particular  interest  since  the  dynamic  loads  in  this  work  are 

applied  normal  to  the  plane  of  the  plate.  The  short  transverse  properties 

for  MIL-S-12560B  armor  steel  have  been  shown  to  differ  significantly 

8 

from  the  properties  in  perpendicular  directions. 

Fractography 

The  IHS  armor  tensile  specimens  fractured  in  a classical  cup-cone 

manner  after  considerable  necking;  the  fracture  appearance  of  the  XAR30 

o 

and  MIL-S-12560B  specimens  was  of  the  45  shear  type,  as  shown  in 
Figure  E-3.  In  the  MIL-S-12560B  steel  large  cracks,  such  as  shown  in 
Figure  E-3(c),  were  often  observed  on  the  specimen  sides.  Figure  E-4  is 
a high  magnification  scanning  electron  micrograph  of  the  central  region 
of  the  cup-cone  of  a specimen  of  improved  homogeneous  steel,  showing 
the  fineness  of  the  dimple  structure.  The  small  spherical  particles 
at  the  centers  of  the  dimples  were  analyzed  by  nondispersive  X-rays  as 
aluminum,  and  are  thus  thought  to  be  aluminum  oxide.  This  fine  uniform 
dispersion  of  these  particles  probably  contributes  to  the  excellent 
strength  and  ductility  characteristics. 

A plate  of  IHS  armor  that  had  been  tested  ballistically  at  AMMRC 
was  sent  to  SRI  for  examination.  We  sectioned  the  plate  through  the 
center  of  the  hole  made  by  the  ballistic  round  and  polished  this  surface 
for  observation  with  an  optical  microscope.  Figure  E-5(a)  shows  the 
fracture  damage  that  occurred  during  projectile  penetration.  It  is 
important  to  note  that  the  cracks  are  rather  randomly  oriented  and  thus 
can  Join  up  to  isolate  fragments,  a process  that  had  indeed  occurred 
at  a location  included  in  the  figure.  Etching  revealed  dense  networks 
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FIGURE  E-4  SCANNING  ELECTRON  MICROGRAPH  OF 

CENTRAL  REGION  OF  CUP-CONE  FRACTURE 
SURFACE  OF  IMPROVED  HOMOGENEOUS 
STEEL  TENSILE  SPECIMEN 
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IMPACT 

DIRECTION 


PROJECTILE 

IMPACT 

DIRECTION 


<•>  POLISHED  SECTION  THROUGH 
PROJECTILE  HOLE 


ETCHED  CLOSE-UP  VIEW  OF  (*) 


MP-2024-6A 


FIGURE  E-5 


PLANAR  MICROCRACKS  AND  ADIABATIC  SHEAR  ZONES 
NEAR  THE  PLUGGED  REGION  IN  IMPROVED  HOMOGENEOUS 
STEEL  ARMOR 


* 

of  adiabatic  shear  bands.  Figure  E-5(b)  is  a closeup  view  of  Figure 
E-5(a),  and  shows  that  many  of  the  cracks  lie  within  the  shear  bands. 
Figure  E-5  demonstrates  the  important  role  of  adiabatic  shear  in  the 
penetration  process  for  IHS,  The  probable  time  sequence  of  events  during 
penetration  is  as  follows: 

(1)  Under  ballistic  impact,  a network  of  adiabatic  shear  bands 
forms  in  the  armor  plate. 

(2)  The  narrow  shear  bands  are  preferred  fracture  sites,  and  planar 
microcracks  nucleate  and  grow  within  them. 

(3)  The  cracks  join  up  with  one  another  and  with  the  surfaces  of 
the  armor  plate,  freeing  the  plug  and  a number  of  fragments. 


Adiabatic  shear  bands  are  narrow  regions  of  highly  localized  large 
plastic  shear  strains.  The  heating  accompanying  the  shear  deformation 
may  increase  the  temperature  to  where  solid  phase  transformations  or 
even  melting  occurs.  Very  rapid  quenching  of  this  material  follows 
because  the  large  volume  of  adjacent  material  in  intimate  contact  with 
the  narrow  band  conducts  the  heat  away  at  high  rates.  In  the  present 
instance  a transformation  to  austenite  followed  by  another  rapid 
transformation  to  martensite  probably  occurred.  The  white  etching 
response  in  5%  nital  is  consistent  with  the  existence  of  martensite. 
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Appendix  F 


THE  MICROMECHANISM  OF  DYNAMIC  FRACTURE  IN  ROLLED  ARMOR  STEEL 

Observations  of  polished  cross  sections  of  shock  damaged  specimens 

of  ferrous  materials  indicate  two  distinctive  fracture  morphologies. 

Rolled  homogeneous  XAR30  armor  steel  for  example,  exhibits  a relatively 

small  number  of  large  parallel  crack-like  artifacts  (as  shown  earlier  in 

3 

Figure  6),  whereas  Armco  iron  acquires  many  more  microcracks,  which  are 
generally  much  shorter,  more  randomly  oriented,  and  more  slit-like,  as 
shown  in  Figure  F-l. 

The  observerations  in  Armco  iron  are  attributable  to  cleavage-type 
fracture  of  individual  grains.  The  grains  and  hence  the  cleavage  planes 
are  nearly  randomly  oriented,  which  explains  the  observed  random  micro- 
crack orientation  and  also  the  limited  size  of  the  microcracks — the  latter 
because  large  deviations  from  crystallographic  match-up  between  grains 
makes  it  difficult  for  a cleavage  crack  to  form  in  an  adjacent  grain. 

Thus  most  microcracks  in  Armco  iron  arrest  at  grain  boundaries. 

But  whereas  the  fracture  morphology  in  Armco  iron  is  controlled  by 
the  presence  of  weak  cleavage  planes,  the  morphology  of  rolled  steels 
reflects  the  weak  planes  of  the  rolling  direction.  Examination  of  the 
large  planar  fractures  in  XAR30  armor  steel  at  high  magnification  reveals 
that  each  large  fracture  comprises  a large  number  of  microfractures 
having  a spherical  appearance.  It  appears  that  the  microfractures  nucleate 
in  a plane,  grow  spherically,  and  coalesce  with  one  another  to  result 
in  the  observed  planar  macrofractures. 

Full  spall  or  the  formation  of  a continuous  macrofracture  within  a 
specimen  occurs  when  adjacent  planar  fractures  coalesce.  The  micrograph 
in  Figure  F-2  shows  two  parallel  but  nonplanar  fractures  in  the  process 
of  coalescing.  A profusion  of  tiny  microfractures  has  formed  in  the 
path  of  incipient  coalescence  under  the  action  of  high  shear  stresses, 
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FIGURE  F-1  POLISHED  CROSS  SECTION  IN  SHOCK 
LOADED  ARMCO  IRON  REVEALING 
INTERNAL  CLEAVAGE  CRACKS 
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FIGURE  F-2  POLISHED  CROSS  SECTION  IN  XAR30  ARMOR  STEEL  SHOWING  INCIPIENT  COALESENCE  OF  TWO 
PLANAR  FRACTURES 


as  in  the  formation  of  the  planar  fractures,  these  microfractures 
nucleate,  grow,  and  join  up  until  a continuous  fracture  connects  the 
two  planes. 

This  process  is  depicted  schematically  in  Figure  F-3. 

The  micrographs  of  shock  loaded  XAR30  armor  steel  in  Figure  F-4 
clearly  show  that  inclusions  are  involved  in  the  nucleation  of  fracture. 
This  leads  us  to  speculate  that  the  improved  ballistic  performance  of 
unidirectionally  solidified  or  electroslag  remelted  steel  is  connected 
with  its  lower  inclusion  content. 

The  view  in  Figure  F-4  is  slightly  etched  to  reveal  the  grey  sulfide 
inclusions  but  not  the  grain  structure.  The  inclusions,  whose  elongated 
form  and  biased  orientation  result  from  the  rolling  process,  have 
fractured  in  a brittle  manner.  As  these  cracks  attempt  to  extend  into 
the  more  ductile  steel,  their  morphology  becomes  more  equiaxed.  Figure 
F-4(b)  is  a more  heavily  etched  area  of  the  specimen,  showing  a large 
planar  fracture  that  apparently  originated  at  the  broken  inclusion, 
extended  by  coalescence  of  more  ductile  microfractures,  and  finally 
coalesced  with  a parallel  but  nonplanar  fracture  by  a similar  process. 

The  apparent  nests  of  voids  observable  at  high  magnifications  may 

not  be  caused  by  the  same  dislocation  mechanisms  responsible  for  void 

formation  in  aluminum,  copper,  and  tantalum,  such  as  described  by  Stevens 
18,19 

et  al.  Rather,  as  implied  from  dynamic  fracture  work  in  Armco  iron, 

the  "voids"  may  nucleate  as  planar  microcracks  on  cleavage  planes,  which 
upon  reaching  the  grain  boundaires,  arrest  and  begin  to  widen  (Figure  F-5) . 
The  result  is  a void-like  microfracture,  which  attained  its  morphology 
by  extensive  plastic  flow  at  a planar  crack  front  instead  of  by  main- 
taining a pseudo-spherical  morphology  during  its  growth  phase. 
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SPHERICAL  EXPANSION  OF  THE  VOIDS  AND  COALESCENCE  OF  THE  ROW.  OTHER  ROWS 
OF  VOIDS  ARE  NUCLEATED  ON  EITHER  SIDE  OF  MAIN  DAMAGED  ZONE. 
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ccco 

COCO  CCCO 


( 3 ) JOINING  UP  OF  THE  ROWS  BY  SHEAR  CRACKS,  WHICH  ALSO  FORM  BY  NUCLEATION.  GROWTH. 
AND  COALESCENCE  OF  MICROVOIDS. 
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FURTHER  WIDENING  AND  COALESCENCE. 


THE  FRACTURE  DAMAGE  OF  STAGE  2 APPEARS  MACROSCOPICALLY  TO  BE  PLANAR  MICROCRACKS 
FRACTURE  MODE  IN  ACTUALITY  IS  A DUCTILE,  ENERGY-ABSORBING  MODE.  SEM  WORK  SHOWS 
NESTS  OF  VOIDS  CONNECTED  BY  SHEAR  CRACKS  LYING  AT  STEEP  ANGLES. 


MANY  INCLUSIONS;  MANY  HAVE  CRACKS  PRODUCED  EITHER  DURING  ROLLING  OR  DURING  IMPACT. 
TWO  MAIN  TYPES: 

(a)  GREY  LENSES  (STRINGERS),  PROBABLY  SULFIDE  / J 

(b)  ORANGE  SQUARES  □ 

MA-2024-34 


FIGURE  F-3  PROPOSED  DYNAMIC  FRACTURE  PROCESS  IN  ROLLED  ARMOR  STEEL 
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(a)  POLISHED  SPECIMEN,  1000X 


(b)  ETCHEO  SPECIMEN,  SOOX 

MP-2024-3B 


FIGURE  F-4  POLISHED  AND  ETCHED  CROSS  SECTION 
OF  SHOCK  LOADED  XAR30  ARMOR  STEEL 
SHOWING  CRACKING  AT  INCLUSIONS 


(i)  CRACK-FREE  GRAIN  BEFORE  IMPACT 


(c)  ARRESTED  CRACK  BEGINS  TO  WIOtN 
AS  THE  MATERIAL  AT  THE  CRACK  TIP 
UNDERGOES  INTENSE  PLASTIC  DEFORMATION 


(b)  PLANAR  CLEAVAGE  CRACK  IS 

PRODUCED  DURING  IMPACT  AND 
ARRESTS  AT  GRAIN  BOUNDARY 


Id)  WIDENING  HAS  OCCURRED  TO  SUCK 
AN  EXTENT  THAT  THE  ORIGINAL 
PLANAR  FRACTURE  APPEARS  VOID-LIKE 

MP-2024-36 


FIGURE  F-5  PROPOSED  MECHANISM  FOR  VOID  FORMATION  IN  ARMOR 
STEEL  UNDER  DYNAMIC  LOADING 
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The  suggested  mechanism  for  void  formation  in  armor  steel  is 
nucleation  in  the  form  of  planar  microcracks,  arrest  at  grain  boundaries, 
and  displacement  of  crack  faces.  In  coarse  grained  material  such  as 
Armco  iron,  the  planar  crack  morphology  is  still  apparent  even  after 
significant  widening,  as  shown  in  Figure  F-5.  However,  for  steels  in 
which  the  average  grain  diameter  is  about  the  same  as  the  crack  tip 
stret'’1'  the  planar  fractures  of  the  original  microcrack  are  lost  and 
the  ,..^cr'>f racture  appears  void-like. 

The  dynamic  fracture  process  in  rolled  armor  steel  is  based  on  the 
proposed  void  formation  mechanism  and  is  thought  to  occur  in  the  manner 
indicated  in  Figure  F-3.  Large  numbers  of  inclusions  fracture  and/or 
voids  nucleate  on  planes  parallel  to  the  rolling  plane.  These  fractures 
become  more  equiaxed  as  they  grow,  and  coalesce  to  form  large  planar 
fractures.  Shear  cracks  form  between  parallel  nonplanar  macrocracks. 
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